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ABSTRACT 

We present a multiwavelength analysis of a sample of four hot (T x > 8keV) X-ray galaxy clusters (A1689, 
A2261, A2142, and A2390) using joint AMiBA Sunyaev-Zel'dovich effect (SZE) and Subaru weak lensing 
observations, combined with published X-ray temperatures, to examine the distribution of mass and the intr- 
acluster medium (ICM) in massive cluster environments. Our observations show that A2261 is very similar 
to A1689 in terms of lensing properties. Many tangential arcs are visible around A2261, with an effective 
Einstein radius ~ 40" (at z ~ 1 .5), which when combined with our weak lensing measurements implies a mass 
profile well fitted by an NFW model with a high concentration c v ; r ~ 10, similar to A 1689 and to other massive 
clusters. The cluster A2142 shows complex mass substructure, and displays a shallower profile (c v j r ~ 5), con- 
sistent with detailed X-ray observations which imply recent interaction. The AMiBA map of A2142 exhibits 
an SZE feature associated with mass substructure lying ahead of the sharp north-west edge of the X-ray core 
suggesting a pressure increase in the ICM. For A2390 we obtain highly elliptical mass and ICM distributions 
at all radii, consistent with other X-ray and strong lensing work. Our cluster gas fraction measurements, free 
from the hydrostatic equilibrium assumption, are overall in good agreement with published X-ray and SZE 
observations, with the sample-averaged gas fraction of (/ gas (< rzoo)) = 0.133 ±0.027, for our sample with 
(M vir ) = (1.2±0.1) x IO^Mq/T 1 . When compared to the cosmic baryon fraction = Qj/ Q, m constrained by 
the WMAP 5-year data, this indicates (/ ga s,20o)//i> = 0.78 ±0.16, i.e., (22 ± 16)% of the baryons are missing 
from the hot phase of clusters. 

Subject headings: cosmology: observations, cosmic microwave background — galaxies: clusters: individual 
(A1689, A2142, A2261, A2390) — gravitational lensing 



1. INTRODUCTION 

Clusters of galaxies, the largest virialized systems known, 
are key tracers of the matter distribution in the large scale 
structure of the Universe. In the standard picture of cosmic 
structure formation, clusters are mostly composed of dark 
matter (DM) as indicated by a great deal of observational evi- 
dence, with the added assumptions that DM is non relativistic 
(cold) and collisionless, referred to as CDM. Strong evidence 
for substantial DM in clusters c omes from multiwavele ngth 
studies of interacting clusters dMarkevitch et alJ 120021) . in 
which weak gravitational lensing of background galaxies en- 
ables us to directly map the distribution of gravitating matter 
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in merging clust ers regardless of the physical/dynamical stat e 
of the system dClowe et al.1 l2006t lOkabe & Umetsu! 120081) . 
The bulk of the baryons in clusters, on the other hand, reside 
in the X-ray emitting intracluster medium (ICM), where the 
X-ray surface brightness traces the gravitational mass domi- 
nated by DM. The remaining baryons ar e in the form of lumi- 
nous galaxies and fain t intracluster light (Fukugita et al.ll998t 
Gonza lez et al.l l2005). Since rich clusters represent high den- 
sity peaks in the primordial fluctuation field, their baryonic 
mass fraction and its redshift dependence can in p rinciple 
be used to constrain the background cosmology (e.g., Sasaki 
1996t lAJlen etalJl200l 12004 120081) . In particular, the gas 
mass to total mass ratio (the gas fraction) in clusters can 
be used to place a lower limit on the cluster baryon frac- 
tion, which is expected to match the cosmic baryon fraction, 
fi, = f2j/fi m . However, non-gravitational processes associ- 
ated with cluster formation, such as radiative gas cooling and 
AGN feedback, would break the self-similarities in cluster 
properties, which can cause the gas fraction to acquire so me 
mass dependence dBialek et al.ll200UlKravtsov et alj|2005l) . 

The deep gravitational potential wells of massive clus- 
ters generate weak shape distortions of the images of back- 
ground sources due to differential deflection of light rays, 
resulting in a systematic distortion pattern around the cen- 
ters of massive clusters, known as weak gravitational lens - 
ing (e.g. JUmetsu et al.ll999tlBartermann & Schneiderl2001l) . 
In the past decade, weak lensing has become a power- 
ful, reliable measure to map the distribution of matter in 
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clusters, dominated by invisible DM, without requiring any 
assumption about the physical/dynamical state of the sys- 
tem (e.g., lOoweitiD EOQl |Ok*iAUmiti3l2008). Re- 
cently, cluster weak lensing has b een used to examine the 
form of DM density profiles (e.g., iBro adhurst et al. 2005bl 
I2008t iMandelbaum et al]l2008t lUmetsu & Broadhursdl2008h . 
aiming for an observational test of the equilibrium density 
profile of DM halos and the scaling relation between halo 
mass and concentration, predicted by A^-body simulations 
in the standard Lambda Cold Dark Mat ter (ACDM) model 
( Sperg eTet alj|2007t [Komatsu et al. 2008). Observational re- 
sults show that the form of lensing profiles in relaxed clusters 
is consistent with a continuously steepening density profile 
with in creasing radius, wel l described by the general NFW 
model (Navarro et alj|1997h . expected for collisionless CDM 
halos. 

The Yuan -Tseh Lee Arr ay for Microwave Background 
Anisotropy dHo et alj|2008l) is a platform-mounted interfer- 
ometer array of up to 19 elements operating at 3 mm wave- 
length, specifically designed to study the structure of the cos- 
mic microwave background (CMB) radiation. In the course of 
early AMiBA operations we conducted Sunyaev-Zel'dovich 
effect (SZE) observations at 94 GHz towards six massive 
Abell clusters with the 7 -element compact array dWu et al.l 
l2008al) . At 94 GHz, the SZE signal is a temperature decre- 
ment in the CMB sky, and is a measure of the thermal 
gas pressure in the ICM integrat ed along the line of sight 
dBirkinshawll 19991; iRephaelill 19951) . Therefore it is rather in- 
sensitive to the cluster core as compared with the X-ray data, 
allowing us to trace the distribution of the ICM out to large 
radii. 

This paper presents a multiwavelength analysis of four 
nearby massive clusters in the AMiBA sample, A 1689, 
A2261, A2142, and A2390, for which high-quality deep 
Subaru images are available for accurate weak lensing mea- 
surements. This AMiBA lensing sample represents a sub- 
set of the high-mass clusters that c an be selected by their 
high (T x > 8keV) gas temperatures dWu et all 12008 a). Our 
joint weak lensing and SZE observations, combined with sup- 
porting X-ray information available in the published litera- 
ture, will allow us to constrain the cluster gas fra ctions with- 
out the assumption of hy drostatic equilibrium (Mye rs et alj 
ll997tlUmetsu et aTll2005l) . complementing X-ray based stud- 
ies. Our companion papers complement details of the instru- 
ments, system performance and verification, observati ons and 
data an alysis, and early science results from AMiBA. lHo et alj 
(2008) describe the design concepts and specifications of the 
AMiBA tele scope. Technical aspe cts of the instrume nts are 
described in lChen et ail d2008) and iKoch et ail d2008al) . De- 
tails of th e first SZE observ ati ons and data analysis are pre- 
sented in IWu et all d2008al) . iNishioka et afl ([2008) assess 
the integ r ity of AMiBA data with several statistical tests. 
iLin et al] d2008l) d i scuss the system performance and verifi- 
cation. " iLiu et alJ d2008l) examine the levels of contamina- 
tion from foregrou nd sources and the primary CMB radiation. 
IKoch et ail (2008b) present a measurement of the Hubble con- 
stant, H Q , from AMiBA SZE and X-ray data. iHuang et al] 
(2008) discuss cluster scaling relations between AMiBA SZE 
and X-ray observations. 

The paper is organized as follows. We briefly summarize in 
§2 the basis of cluster SZE and weak lensing. In §3 we present 
a concise summary of the AMiBA target clusters and observa- 
tions. In §4 we describe our weak lensing analysis of Subaru 
imaging data, and derive lensing distortion and mass profiles 



for individual clusters. In §5 we examine and compare cluster 
ellipticity and orientation profiles on mass and ICM structure 
in the Subaru weak lensing and AMiBA SZE observations. 
In §6 we present our cluster models and method for measur- 
ing cluster gas fraction profiles from joint weak-lensing and 
SZE observations, combined with published X-ray tempera- 
ture measurements; we then derive cluster gas fraction pro- 
files, and constrain the sample-averaged gas fraction profile 
for our massive AMIBA-lensing clusters. Finally, a discus- 
sion and summary are given in §7. 

Throughout this paper, we adopt a concordance 
ACDM cosmology with Q m = 0.3, Q\ = 0.7, and 
h = H()/(100kms"'Mpc" 1 ) = 0.7. Cluster properties are 
determined at the virial radius r vu - and radii (r2oo,'"500,'"250o), 
corresponding to overdensities (200, 500, 2500) relative to the 
critical density of the universe at the cluster redshift. 

2. BASIS OF CLUSTER SUNYAEV-ZEL'DOVICH EFFECT AND 
WEAK LENSING 

2.1. Sunyaev-Zel 'dovich Effect 

We begin with a brief summary of the basic equations of the 
thermal SZ E. Our notation here closely follows the standard 
notation of Re phaelil d 19951) . 

The SZE is a spectral distortion of the CMB radiation re- 
sulting from the inverse Compton scattering of cool CMB 
photons by the hot ICM. The non-relativistic form of the spec- 
tral change was obtained by Sunyaev-Zel'dovich (1972) from 
the Kompaneets equation in the non-relativistic limit. The 
change in the CMB intensity /cmb due to the SZE is writ- 
ten in terms of its spectral function g and of t he integral of the 
electron pressure along the line-of-sigh t as ( Rephaelil 11995b 
lBirkinshaw|[T999HCarlstrom et alj|2002l) : 



A/ SZE (i/) = *norm 



(1) 



where x(v) is the dimensionless frequency, x = 
hu/(k B T CM B) ~ 1.66(^/94 GHz), with k B being 
the Boltzmann constant and 7cmb = 2.725 K be- 
ing the CMB temperature at the present-day epoch, 
/no,™ = (2h / c 3 ) (k B T CMB / h) 2 ~ 2.7 x lO^ysr" 1 , and y(8) is 
the Comptonization parameter defined as 



y = 



dl (T t h n e 



k R T„ 



m e c- 



dl — — ksT„ as , 
fj, e m p 



(2) 

where er,/,, m e , c, and \i e are the Thomson cross section, the 
electron mass, the speed of light, and the mass per electron 
in units of proton mass m p , respectively; for a fully ionized 
H-He plasma, fj, e = 2/(1 +X) ~ 1.14, with X being the Hydro- 
gen primordial abundance by mass fraction, X ~ 0.76; r max is 
the cutoff radius for an isolated cluster (see 36.31 ). The SZE 
spectral function g(x) is expressed as 



g(x) = gNR(x) [l + $sze(x, T gl 



(3) 



where ^nrM is the thermal spectral function in the non- 
relativistic limit (Sunyaev & Zel'dovich 1972), 



gm(x) = 



4 r 

x e 



(e*-l) 2 V e x -\ 



e x +\ „ 
x --4 



(4) 



which is zero at the cross-over frequency xq ~ 3.83, or 
vq = 217 GHz, and 5§ze(x, Teas) is the relati vistic correction 
dChallinor & Lasenbyl fl998; Itoh etlii] [19981) . The fractional 
intensity decrease due to the SZE with respect to the pri- 
mary CMB is maximized at v ~ 100 GHz (see Figure 1 of 



Umetsu et al. 



3 



Zhan g et alj|2002lh which is well matched to the observing 
frequency range 86-102 GHz of AMiBA. At the central fre- 
quency v c = 94 GHz of AMIBA, g(x) ~ -3.4. For our hot X- 
ray clusters with Tx = 8-10keV, the relativistic correction to 
the thermal SZE is 6-7% at v c = 94 GHz. 

2.2. Cluster Weak Lensing 

Weak gravitational lensing is responsible for the weak 
shape-distortion and magnification of the images of back- 
ground sources due to the gravitational field of intervening 
foreground clusters of galaxies and large scale structures in 
the universe. The deformation of the image can be described 
by the 2x2 Jacobian matrix A a p (a, (3 = 1 , 2) of the lens map- 
ping. The Jacobian A a p is real and symmetric, so that it can 
be decomposed as 

A a f3 = (l-K)5 al 3-r al3 , (5) 



_/+7i 72 
Taf3 -{ 72-7i 



(6) 



where 6 a p is Kronecker's delta, r Q ^ is the trace-free, sym- 
metric shear matrix with j a being the components of spin- 
2 complex gravitational shear 7 := 71 + £72, describing the 
anisotropic shape distortion, and k is the lensing convergence 
responsible for the trace-part of the Jacobian matrix, describ- 
ing the isotropic area distortion. In the weak lensing limit 
where k, \ j\ <C 1, T a/ 3 induces a quadrupole anisotropy of the 
background image, which can be observed from ellipticities 
of background galaxy images. The flux magnification due to 
gravitational lensing is given by the inverse Jacobian determi- 
nant, 

1 1 m 

M deU (1- K ) 2 -| 7 | 2 ' ' 

where we assume subcritical lensing, i.e., det4(0) > 0. 

The lensing convergence is expressed as a line-of-sight pro- 
jection of the matter density contrast S m = (p m — p)/p out to 
the source plane (s) weighted by certain combination of c o- 
moving angular diameter distances r (e.g. JJain et alJ uOOO). 

5, 



K = 



3//q £l m 



2c 2 jo 

KxMxs-x) 



dxG(X,Xs)—= ldZ m Z-±, (8) 



r(Xs) 



(9) 



where a is the cosmic scale factor, x i s me co-moving 
distance, £,„ is the surface mass density of matter, E m = 
J dx<*(pm~ P)i with respect to the cosmic mean density p, and 
Edit is the critical surface mass density for gravitational lens- 



ing, 



4ttG D d D ds 



(10) 



with D s , Dj, and Dj, being the (proper) angular diameter dis- 
tances from the observer to the source, from the observer to 
the deflecting lens, and from the lens to the source, respec- 
tively. For a fixed background cosmology and a lens redshift 
Zd, E cl it is a function of background source redshift z s - For 
a given mass distribution £ m (0), the lensing signal is propor- 
tional to the angular diameter distance ratio, D c / S /D s . 

In the present weak lensing study we aim to reconstruct 
the dimensionless surface mass density k from weak lens- 
ing distortion and magnification data. To do t his, we utilize 
the relation between t he gradients of k and 7 ( Kaiser] [l 9951: 
ICrittenden et~aT1l2002l) . 

An(9) = d a d< 3 T af3 (9) = 2£>* 7 (0) (11) 



where t> is the complex differential operator t> = (d 2 -d 2 )/2+ 
id\d2- The Green's function for the two-dimensional Pois- 
son equation is A~ l (9, 0') = In \9— 0'\/(2tt), so that equation 
( fTTT i can be solved to yield the following non-local relation 
between n and 7 (Kaiser & Squires 1993): 

«(0) = i Jd 2 e'D*(e-e')j(e') (12) 



where D(9) is the complex kernel defined as 
d 2 . — 9 2 — 2.i0\ 62 



D(9): 



\6\< 



(13) 



Similarly, the spin-2 shear field can be expressed in terms of 
the lensing convergence as 



1 



7 (0) = - / d 2 9'D(9-9')n(9'). 



(14) 



Note that adding a constant mass sheet to k in equation ([Pil l 
does not change the shear field 7(0) which is observable in 
the weak lensing limit, leading to the so-called mass-sheet de- 
generacy (see eq. fT6l ) ba sed solely on shape-distortion mea- 
surem ents (e.g., Bartelma nn & Schneiderl200ll : IUmetsu et alJ 
1999). In general, the observable quantity is not the gravita- 
tional shear 7 but the reduced shear, 

(15) 



l-K 



in the subcritical regime where det4 > (or l/g* in the neg- 
ative parity region with det4 < 0). We see that the reduced 
shear g is invariant under the following global transformation: 



k(0)->Ak(0) + 1-A, 7 (0)^A 7 (0) 



(16) 



with an arbitrary scalar constant A y (ISchneider & Seitzl 
fl995h . 

3. AMIBA SUNYAEV-ZEL'DOVICH EFFECT OBSERVATIONS 

3.1. AMIBA Telescope 

The AMiBA is a dual channel 86-102 GHz (3-mm wave- 
length) interferometer array of up to 19-elements with dual 
polarization capabilities sited at 3396 m on Mauna-Loa, 
Hawaii (latitude: +19.5°, longitude: -155.6°) 10 . AMiBA 
is equipped with 4-lag analog, broadband (16 GHz bandwidth 
centered at 94 GHz) correl ators which outpu t a set of 4 real- 
number correlation signals dChen et al.ll20 08). These four de- 
grees of freedom (dof) correspond to two complex visibili- 
ties in two frequency channels. The frequency of AMiBA 
operation was chosen to take advantage of the optimal fre- 
quency window at 3 mm, where the fractional decrement in 
the SZE intensity relative to the primary CMB is close to 
its maximum (see 32.11 ) and contamination by the Galac- 
tic synchrotron emission, dust foregrounds, and the popula- 
tion of cluster/background radio sou rces is minimiz ed (see 
for detailed contamination analysis, iLiu et al.l l2008h . This 
makes AMiBA a unique CMB/SZE interferometer, and also 
complements the wavelength coverage of other existing and 
planned C MB instruments: i nterferometers su ch as AMI 
at 15 GHz dKneissl et al.ll2001l) CBI at 30GHz fPadin et alJ 
I2001L I2002t iMason et all 120031: iPearson et all 120031). SZA 
at 30 and 90 GHz dMroczkowskiet al.ll2008l) . and VSA 11 at 

10 http://amiba.asiaa.sinica.edu.tw/ 

1 1 http://astro.uchicago.edu/sza/ 
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34 GHz ( Watso n et al.l2003l): bolom eter arrays such as ACT, 12 
APEX-SZ^ lHalverson et ani2008l) . and SPT. 14 

In the initial operation of AMiBA, we used seven 0.6 m 
(0.58 m to be precise) Cassegrain antennas dKoch et al.ll2006l) 
co-mounted on a 6 m hex apod platform i n a hexagonal close- 
packed configuration (see lHo et al.ll2008l) . At each of the fre- 
quency channels centered at about 90 and 98 GHz, this com- 
pact configuration provides 21 simultaneous baselines with 
three baseline lengths of d = 0.61, 1.05, and 1.21m, corre- 
sponding to angular multipoles I = 2irV u 2 + v 2 (= 2ird/\) of 
/ « 1 194,2073,2394 at v c = 94GHz. This compact 7-element 
array is sensitive to multipole range 800 / 2600. With 
0.6-m antennas, the instantaneous fie ld-of-view of AMiBA is 
about 23' FWHM (IWu et al.ll2008al) . and its angular resolu- 
tion ranges from 2' to 6' depending on the configuration and 
weighting scheme. In the compact configuration, the angu- 
lar resolution of AMiBA is about 6' FWHM using natural 
weighting (i.e., inverse noise variance weigh ting). The poin t 
source sensitivity is estimated to be ~ 63 mJy dLin et al. 2008) 
in 1 hour of on-source integration in 2-patch main-trail/lead 
differencing observations, where the overall noise level is in- 
creased by a factor of \fl due to the differencing. 

3.2. Initial Target Clusters 

The AMiBA lensing sample, A1689, A2142 A2261 
A2390 , is a subset of the AMiBA cluster sample (see lWu et alj 
2008a), composed of four massive clusters at relatively low 
redshifts of 0.09 £ z ~ 0.23 with the median redshift of z ~ 
0.2. The sample size is simply limited by the availability of 
high quality Subaru weak lensing data. A 1689 is a relaxed, 
round system, and is one of the best studied clusters for lens- 
ing work (e.g.,[B roadhurst et al. 2005b; Limousin et al.ll2007l : 
lUmetsu & Broadhursdl2"00"8t IBroadhurst et al.ll2008l) . A2261 
is a compact cluster with a regular X-ray morphology. A2142 
is a merging cluster with two sharp X-ray surface brightness 
discontinuities in the cluster core (Mar kevitch etail [2000; 
lOkabe & Umetsul 120081) . A2390 shows an elongated mor- 
phology both in the X-ray emission and strong-lensing mass 
distributions dAllen et aU200UlFrve & Broadhurstl 1998b . Ta- 
ble Q] summarizes the physical properties of the four target 
clusters in this multiwavelength study. 

In 2007, AMiBA with the seven small antennas (hence- 
forth AMiBA7) was in the science verification phase. For 
our initial observations, we therefore selected those tar- 
get clusters observable from Mauna Loa during the observ- 
ing period that were known to have strong SZEs at rel- 
atively low redshifts (0.1 ^z^,0.3) from previ ous experi- 
ments , such as OVRO observations at 30 GHz (Mas on et al. | 
1 200 ll) BIMA/OVRO obs ervations at 30 GHz (Gre go etal.1 
I2OO 1 at iReeseetalJ l2002h . VSA observations at 34 GHz 
dLancaster et all2005l). and SuZIE II observations at 145, 221, 
and 355 GHz dBenson et al.l2004l) . The targeted redshift range 
allows the target clusters to be resolved by the 6' resolution 
of AMiBA7, allowing us to derive useful measurements of 
cluster SZE profiles for our multiwavelength studies. At red- 
shifts of z ^ 0.3 (0.2), the angular resolution of AMiBA7 cor- 
responds to £j 560kpc/r' (~ 400kpc/i _1 ) in radius, which is 
< 30^10% (~ 20-30%) of the virial radius 1.5-2MPC/T 1 of 
massive clusters. The requirement of being SZE strong is to 
ensure reliable SZE measurements at 3 mm with AMiBA7. 

12 http://www.hep.upenn.edu/act/act.html 

1 3 http ://bolo . berkeley. edu/apexsz 

14 http://pole.uchicago.edu 



We note that AMiBA and SZA are the only SZE instru- 
ments measuring at 3 mm, but complimentary in their baseline 
coverage. With sensitivities of 20-30 mJy beam" 1 typically 
achieved in 2-patch differen cing observations in 5-10 hours 
of net on-source integration (IWu et al.ll2008al) . we would ex- 
pect ^ 5a detections of SZE fluxes ^ 100-150mJy at 3 mm. 
Finally, our observing period (April-August 2007) limited 
the range of right ascension (RA) of targets, 15 since we re- 
stricted our science observations to nights (roughly 8pm to 
8am), where we would expect high gain stabil ity because the 
ambie nt temperature varies slowly and little (Nishio ka et al.1 
2008). The SZE strong clusters in our AMiBA sample are 
likely to have exceedingly deep potential wells, and indeed 
our AMiBA sample represents a class of hot X-ray clusters 
with observed X-ray temperatures exceeding 8keV (see Ta- 
ble Q]). We note that this may affect the generality of the 
results presented in this study. A main-trail/lead differenc- 
ing scheme has been used in our targeted cluster observa- 
tions where the trail/lead (blank) field is subtracted from the 
main (cluster) field. This differencing scheme sufficiently re- 
moves contamination from groun d spillover and e lectronic 
DC offset in the correlator output dWu et alj|2008al) . A full 
description of AMiBA observations and analysis of the initial 
six target clusters, including the observation strategy, analysis 
methodology, c alibrations, and map-making, can be found in 
IWuet alJd2008allbh . 

4. SUBARU WEAK LENSING DATA AND ANALYSIS 

In this section we present a technical description of our 
weak lensing distortion analysis of the AMiBA lensing sam- 
ple based on Subaru data. The present work on A 1689 is 
based o n the same Subaru image s as analyzed in our earlier 
work of IBroadhurst et alJ d2005bl) and lUmetsu & Broadhursl 
(2008), but our improved color selection of the red back- 
ground has increased the sample size by ~ 16% Q4.3t . This 
work on A2142 i s base d on the same Subaru images as in 
Oka be & Ume tsu (2008), but our inclusion of blue, as well 
as red, galaxies has increased the sample size by a factor 
of 4 (cf. Table 6 of Okabe & Umetsu 2008), leading to a 
significant improvement of our lensing measurements. For 
A2261 and A2390 we present our new weak lensing analysis 
based on Suprime-Cam imaging data retrieved from the Sub- 
aru archive, SMOKA. The read er only interested in the main 
result may skip directly to jj4.4l 

4.1. Subaru Data and Photometry 

We analyze deep images of four high mass clusters in the 
AMiBA sa mple taken by the wide -field camera Suprime-Cam 
(34' x 27'; iMivazaki etaT] 120021) at the prime-focus of the 
8.3m Subaru telescope. The clusters were observed deeply 
in two optical passbands each with seeing in the co-added 
images ranging from 0.55" to 0.88" (see Table |2|. For 
each cluster we select an optimal combination of two filters 
that allows for an efficient separation of cluster/background 
galaxies based on color-magnitude correlations (see Table 
|2]). We use either R c or i' band for our weak lensing mea- 
surements (described in 34.2b for which the instrumental re- 
sponse, sky background and seeing conspire to provide the 
best-quality images. Th e standard pipeli ne reduction soft- 
ware for Suprime-Cam ( Yag i et all 120021) is used for flat- 
fielding, instrumental distortion correction, differential refrac- 
tion, sky subtraction and stacking. Photometric catalogs are 

15 The elevation limit of AMiBA is 30°. 
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con structed from stacked an d matched images using SExtrac- 
tor dBertin & Arnoutsll 19961) . and used for our color selection 
of background galaxies (see 34.31 ). 



4.2. Weak Lensing Distortion Analysis 

We use the IMCAT package developed by N. Kaiser 16 
to perform object detection, photometry and shape mea- 
surem ents, following the formalism outlined in Kais er et"aT] 
(1 19951 KSB). Our analysis pi peline is impleme nted based on 
the procedures described in lErben et al.1 (1200 ll) and on verifi- 
cation tests with STEP1 da ta of mock ground-based observa- 
tions dHevmans et al. 2006). The same analysis pipeline has 
been used in Umetsu & Broadhurst (2008), Okabe & Umetsu 
d2008l) . and lBroadhurst et al.1 d2008l) . 

4.2.1. Object Detection 

Objects are first detected as local peaks in the image by 
using the IMCAT hierarchical peak-finding algorithm hfind- 
peaks which for each object yields object parameters such as 
a peak position ( x), an estimate of the object size (r g ), the sig- 
nificance of the peak detection {v). The local sky level and its 
gradient are measured around each object from the mode of 
pixel values on a circular annulus defined by inner and outer 
radii of 16 x r g and 32 x r g . In order to avoid contamina- 
tion in the background estimation by bright neighboring stars 
and/or foreground galaxies, all pixels within 3 x r g of another 
object are excluded from the mode calculation. Total fluxes 
and half-light radii (r/,) are then measured on sky-subtracted 
images using a circular aperture of radius 3 x r g from the ob- 
ject center. Any pixels within 2.5 x r g of another object are 
excluded from the aperture. The aperture magnitude is then 
calculated from the measured total flux and a zero-point mag- 
nitude. Any objects with positional differences between the 
peak location and the weighted-centroid greater than \d\ = 0.4 
pixels are excluded from the catalog. 

Finally, bad objects such as spikes, saturated stars, and 
noisy detections need to be removed from the weak lensing 
analysis. We removed from our detection catalog extremely 
small or large objects with r g < 1 or r g > 10 pixels, object s 
with low detection significance, v < 7 (see lErben et al.|2001) , 
objects with large raw ellipticities, \e\ > 0.5 (see 34.2.21 1. 
noisy detections with unphysical negative fluxes, and objects 
containing more than 10 bad pixels, nbad > 10. 

4.2.2. Weak Lensing Distortion Measurements 

To obtain an estimate of the reduced shear, g a = 7 a /(l - k) 
(a = 1 , 2), we measure using the getshapes routine in IMCAT 
the image ellipticity e a = {Qu-Qi2,Q\i\ /(Qn + Qii) from 
the weighted quadrupole moments of the surface brightness 
of individual galaxies defined in the above catalog, 



Q a p= d'OWWJem (a,0= 1,2) (17) 

where 1(6) is the surface brightness distribution of an object, 
W(8) is a Gaussian window function matched to the size of 
the object, and the object center is chosen as the coordinate 
origin. In equation (Tl7] > the maximum radius of integration is 
chosen to be # max = 4^. 

Firstly the PSF anisotropy needs to be corrected using the 
star images as references: 



■ P - p a P n * 



' http://www.ifa.hawaii.edu/kaiser/imcat 



(18) 



where P sm is the smear polarizability tensor (which is close 
to diagonal), and q* a = (.P* m )^ge* is the stellar anisotropy ker- 
nel. We select bright, unsaturated foreground stars identified 
in a branch of the half-light radius vs. magnitude diagram to 
measure q* a . In order to obtain a smooth map of q* a which 
is used in equation (fTSI l. we divided the co-added mosaic 
image (of ~ 10K x 8K pixels) into rectangular blocks. The 
block length is based on the coherent scale of PSF anisotropy 
patterns, and is typically 2Kpixels. In this way the PSF 
anisotropy in individual blocks can be well described by fairly 
low-order polynomials. We then fitted the q* in each block 
independently with second-order bi-polynomials, q"{6), in 
conjunction with iterative outlier rejection on each compo- 
nent of the residual: 6e* a = e* a -(P; m ) al3 q* p (6). The final stel- 
lar sample contains typically 500-1200 stars. Uncorrected 
ellipticity components of stellar objects have on average a 
mean offset (from a value of zero) of 1-2% with a few % 
of rms, or variation of PSF across the data field (s ee, e.g., 
lUmetsu & Broadhurstll2008l:IOkabe & Umetsull2008h . On the 
other hand, the mean residual stellar ellipticity <5e* after cor- 
rection is less than or about 10~ 4 , with the standard error on 
this measurement, a few xlO -4 . We show in Figure [T] the 
quadrupole PSF anisotropy fields as measured from stellar el- 
lipticities before and after the anisotropic PSF correction for 
our target clusters. Figure [2] shows the distributions of stel- 
lar ellipticity components before and after the PSF anisotropy 
correction. In addition, we adopt a conservative magnitude 
limit m < 25.5-26.0 ABmag, depending on the depth of the 
data for each clus ter, to avoid systematic err ors in the shape 
measurement (see Umetsu & Broadhurst 2008). From the rest 
of the object catalog, we select objects with r/, > r* h +a(rf t ) 
pixels as a magnitude-selected weak lensing galaxy sample, 
where r* h is the median value of stellar half-light radii r* h , cor- 
responding to half the median width of the circularized PSF 
over the data field, and o-(r* h ) is the rms dispersion of r* h . 

Second, we need to correct image ellipticities for the 
isotropic smearing effect caused by atmospheric seeing and 
the window function used for the shape measurements. The 
pre-seeing reduced shear g a can be estimated from 

g a = (P g X ) af }e' p (19) 

with the pre-seeing shear polarizability tensor P s afi defined as 
iHoekstraetafld 19981) . 



p s 

r afi 



. psh ^psm^psm*-y\psh*~\ 



a/3 ' 



a/3 



nsm 



tr[P 



:sll * 



(20) 



tr[ p sm *] 

with P sh being the shear polarizability tensor; In the sec- 
ond equality we have used a trace approximation to the stel- 
lar shape tensors, P sh * and p sm * . To apply equation ( [T9l i 
the quantity tx[P sh *]/tx[P sm *] must be known for each of the 
galaxi es with different sizescales. Following Hoekst ra et alj 
(1998), we recompute the stellar shapes P sh * and P sm * in 
a range of filter scales r g spanning that of the galaxy sizes 
(r g = [1, 10] pixels). At each filter scale r g , the median 
(tr[P sh *]/tx[P sm *]) over the stellar sample is calculated, and 
used in equation d20b as an estimate of tr[P sh *]/tx[P sm *]. Fur- 
ther, we adopt a scalar correction scheme, namely 



(21) 

dErben et alJ l200lt iHoekstra et all [19981: 
lUmetsu c^B roadhurstT 20081) . In order to suppress artifi- 
cial effects due to the noisy P* estimated for individual 
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galaxies, we apply filtering to raw P? measurements. We 
compute for each object a median value of P* among 
Af-neighbors in the size and magnitude plane to define the 
object parameter space: firstly, for each object, A^-neighbors 
with raw Pf > are identified in the size (r g ) and magnitude 
plane; the median value of P* is then used as the smoothed 

PJ for the object, (PI), and the variance a g of g = g\ +ig2 is 
calculated using equation (19[ . The dispersion a g is used as 
an rms error of the shear estimate for individual galaxies. We 
take N = 30. Finally, we use the estimator g a = e' a j (P^) for 
the reduced shear. 

4.3. Background Selection 

It is crucial in the weak lensing analysis to make a secure 
selection of background galaxies in order to minimize con- 
tamination by unlensed cluster/foreground galaxies and hence 
to make an accurate determination of the cluster mass pro- 
file; otherwise dilution of the distortion signal arises from 
the inclusion of unlensed galaxies, particul arly at small ra- 
dius where the cluster is relat ively dense dBroadhurst et al.1 
l2005bt iMedezinski et al.ll2007l) . This dilution effect is sim- 
ply to reduce the strength of the lensing signal when aver- 
aged over a local ensemble of galaxies, in proportion to the 
fraction of unlensed galaxies whose orientations are randomly 
distributed, thus diluting the lensing signal relative to the ref- 
erenc e background leve l deriv ed from the background popu- 
lation iMedezhikilLlD (120071) . 

To separate cluster members from the background and 
hence minimize the weak lensing dilution, we follow 
an objective backgroun d sel ection method developed by 
IMedezinski et all (120071) and lUmetsu & Broadhurstl (120081) . 
We select red galaxies with colors redder than the color- 
magnitude sequence of cluster E/SO galaxies. The sequence 
forms a well defined line in object color-magnitude space 
due to the richness and relatively low redshifts of our clus- 
ters. These red galaxies are expected to lie in the background 
by virtue of ^-corrections which are greater than for the red 
cluster sequence galaxies; T his has been convincing ly demon- 
strated spectroscopically by iRines & Gelled 12001 . We also 
include blue galaxies falling far from the cluster sequence to 
minimize cluster contamination. 

Figure[3]shows for each cluster the mean distortion strength 
averaged over a wide radial range of 9 = [1', 18'] as a func- 
tion of color limit, done separately for the blue (left) and red 
{right) samples, where the color boundaries for the present 
analysis are indicated by vertical dashed lines for respec- 
tive color samples. Here we do not apply area weighting 
to enhance the effect of d ilution in the central region (see 
lUmetsu &"B roadhurst 2008). A sharp drop in the lensing sig- 
nal is seen when the cluster red sequence starts to contribute 
significantly, thereby reducing the mean lensing signal. Note 
that the background populations do not need to be complete in 
any sense but should simply be well defined and contain only 
background. For A1689, the weak lensing signal in the blue 
sample is systematically lower than that of the red sample, 
so that blue galaxies in A 1689 are excluded from the present 
analysis, as was done in lUmetsu & Broadhurstl (12008): on the 
other hand, our improved color selection for the red sample 
has led to a ~ 16% increase of red galaxies. In the present 
stu dy we use for A2142 the same Subaru images as analyzed 
by lOkabe & Umetsul (12008b . but we have improved signifi- 
cantly our lensing measurements by including blue, as well as 
red, galaxies, where the sample size has been increased by a 
factor of 4. 



An estimate of the background depth is required when con- 
verting the observed lensing signal into physical mass units, 
because the lensing signal depends on the source redshifts in 
proportion to Dj s /D s . The mean depth is sufficient for our 
purposes as the variation of the lens distance ratio, Dd s /D s , 
is slow for our sample because the clusters are at relatively 
low redshifts (zd ~ 0.1-0.2) compared to the redshift range 
of the background galaxies. We estimate the mean depth 
(Dds/D s ) of the combined red+blue background galaxies by 
applying our color-magnitude sele ction to Subaru mu lticolor 
photometry of the H DF-N region dCapak et al.l 12004) or the 
COSMOS deep field dCapak et al.ll2007l) . depending on the 
availability of filters. The fractional uncertainty in the mean 
depth (D,{ s /D s ) for the red galaxies is typically ~ 3%, while 
it is about 5% for the blue galaxies. It i s useful to define 
the distance-equivalent source reds hift z s ,d (IMedezins ki et al.1 
120071; lUmetsu & Broadhurstll2008l) defined as 

' Dels \ _ Dth 

D s I ~ D s 



(22) 



= 70 +0 06 95+ - 79 98 +0 24 



1 OO +0 - 25 for 

0.16' i ' uu -0.16 1UI 



We find z s d 

A1689, A2142, A2261, and A2390, respectively. For the 
nearby cluster A2142 at z — 0.09, a precise knowledge of the 
source redshift is not critical at all for lensing work. The mean 
surface number density (n g ) of the combined blue+red sam- 
ple, the blue-to-red fraction of background galaxies (B/R), the 
estimated mean depth (Dd s /D s ), and the effective source red- 
shift Zs.d are listed in Table [2] 

4.4. Weak Lensing Map-Making 

Weak lensing measurements of the gravitational shear field 
can be used to reconstruct the underlying projected mass den- 
sity field. In the present study, we will use the dilution-free, 
color-selected background sample ( jj4.3l ) both for the 2D mass 
reconstruction and the lens profile measurements 17 . 

Firstly, we pixelize distortion data of background galax- 
ies into a regular grid of pixels using a Gaussian w g (0) oc 
exp[-6» 2 /6»}] with 6 f = FWHM/V41n2. Further we incorpo- 
rate in the pixelization a statistical weight u g for an individual 
galaxy, so that the smoothed estimate of the reduced shear 
field at an angular position is written as 



g a {0)-- 



(23) 



where g a j is the reduced shear estimate of the rth galaxy at 
angular position 0,, and u g j is the statistical weight of rth 
galaxy taken as the inverse variance, u g j = 1 / (er 2 i + a 2 ), with 
a g j being the rms error for the shear estimate of ith galaxy 
(see § | 4.2.21l and a 2 being the softening constant variance 
dHamana et alj|2003l) . We choose a = 0.4, which is a typical 
value of the mean rms a g over the background sample. The 
case with a = corresponds to an inverse-variance weighting. 
On the other hand, the limit a 3> a g j yields a uniform weight- 
ing. We have confirmed that our results are insensitive to the 
choice of a (i.e., inverse-variance or uniform weighting) with 
the adopted smoothing parameters. The error variance for the 
smoothed shear g = g\+ ig2 d23l is then given as 



(EiW«,i"*,i) 



(24) 



17 Okabe & Umetsu (2008) used the magnitude-selected galaxy sample 
in their map-making of nearby merging clusters to increase the background 
sample size, while the dilution-free red background sample was used in their 
lensing mass measurements. 
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where w g .; = w g {0-Oi) and we have used (gajgpj) = 
(1 /2)<r% ti 5*g6fj with and 5fj being the Kronecker's delta. 

We then invert the pixelized reduced-shear field d23l to ob- 
tain the lensing convergence field k{6) using equation $Y2\ . In 
the map-making we assume linear shear in the weak-lensing 
limit, that is, g a = 7 a /(l - k) ~ j a . We adopt the Kaiser 
& Squires inversion method (Kaiser & Squi res! Il 9931) . which 
makes use of the 2D Green function in an infinite space (j j2. 2\ . 
In the linear map-making process, the pixelized shear field is 
weighted by the inverse of the variance d24l i. Note that this 
weighting scheme corresponds to using only the diagonal part 
of the noise covariance matrix, N(6j,9j) = (Ag(8j)Ag(0 ;)), 
which is only an approximation of the actual inverse noise 
weighting in the presence of pixel-to-pixel correlation due to 
non-local Gaussian smoothing. In Table [2] we list the rms 
noise level in the reconstructed k{6) field for our sample of 
target clusters. For all of the clusters, the smoothing scale 9/ 
is taken to be Of = 1' (6>fwhm — 1-665'), which is larger than 
the Einstein radius for our background galaxies. Hence our 
weak lensing approximation here is valid in all clusters. 

In Figure |4] we show, for the four clusters, 2D maps of the 
lensing convergence k(9) = S w ( 0)/S cr j t reconstructed from 
the Subaru distortion data ( {34.41 ). each with the correspond- 
ing gravitational shear field overlaid. Here the resolution of 
the k field is ~ 1.665' in FWHM for all of the four clusters. 
The side length of the displayed region is 22', corresponding 
roughly to the instantaneous field-of-view of AMiBA (~ 23' 
in FWHM). In the absence of higher-order effects, weak lens- 
ing only induces curl-free E-mode distortions, responsible for 
tangential shear patterns, while the B-mode lensing signal is 
expected to vanish. For each case, a prominent mass peak is 
visible in the cluster center, around which the lensing distor- 
tion pattern is clearly tangential. 

Also shown in Figure [4] are contours o f the AMiBA flu x 
density due to the thermal SZE obtained by lWu et"aT1 (|2008a). 
The resolution of AMiBA7 is about 6' in FWHM (@. The 
AMiBA map of A 1689 reveals a bright and compact struc- 
ture in the SZE, similar to the compact and round mass 
distribution reconstructed from the Subaru distortion data. 
A2142 shows an extended structure in the SZE elongated 
along the northwest-southeast direction, consistent with the 
direction of elongatio n of the X-ray halo, w ith its general 
cometary appearance (Marke vitch et al.l 120001) . In addition, 
A2142 shows a slight excess in SZE signals located ~ 10' 
northwest of the cluster center, associated with mass sub- 
structure seen in our lensing k map (Figure 0); This slight 
excess SZE appears extended for a couple of synthesized 
beams, although the per-beam sig nificance level is marginal 
(2— 3a). lOkabe & Umetsul (120081) showed that this northwest 
mass substructure is also associated with a slight excess of 
cluster sequence galaxies, lying ~ 5' ahead of the northwest 
edge of the central X-ray gas core. On the other hand, no 
X-ray counterpart to the northwest substructure was found in 
the X- ray images from Chandr a and XMM-Newton observa- 
tions (lOkabe & Umetsul 12008). A2261 shows a filamentary 
mass structure with unknown redshift, ext ending to the west 
of the cluster core (Maughan et al.1 120081) . and likely back- 
ground structures which coincide with redder galaxy concen- 
trations (see 34.5.2l for details). Our AMiBA and Subaru ob- 
servations show a compact structure both in mass and ICM. 
The elliptical mass distribution in A2390 agrees well with the 
shape seen by AMiBA in the thermal SZE, and is also con- 
sistent with other X-ray and strong lensing work. A quantita- 



tive comparison between the AMiBA SZE and Subaru lensing 
maps will be given in $5] 

4.5. Cluster Lensing Profiles 
4.5.1. Lens Distortion 

The spin-2 shape distortion of an object due to gravita- 
tional lensing is described by the complex reduced shear, 
g = g \ +ig2 (see equation fT5l ). which is coordinate depen- 
dent. For a given reference point on the sky, one can instead 
form coordinate-independent quantities, the tangential distor- 
tion g+ and the 45° rotated component, from linear combina- 
tions of the distortion coefficients g\ and gi as 

^+ = -(^iCos20+g 2 sin20), g x = -(g 2 cos2(f>- gi sin 20), 

(25) 

where <fi is the position angle of an object with respect to the 
reference position, and the uncertainty in the g+ and g x mea- 
surement is er+ = a x = cr g / \/2 = a in terms of the rms error 
a g for the complex shear measurement. In practice, the ref- 
erence point is taken to be the cluster center, which is well 
determined by the locations of the brightest cluster galaxies. 
To improve the statistical significance of the distortion mea- 
surement, we calculate the weighted average of g+ and g x , 
and its weighted error, as 



x(#m)) : 



<r+(6m) = (Tx(0m) = 



\ (Ei%i) 



(26) 
(27) 

(28) 



where the index i runs over all of the objects located within the 
rath annulus with a median radius of 9 m , and u gi i is the inverse 
variance weight for /th object, u gt i = 1/(<t^ ; + a 2 ), softened 
with a = 0.4 (see ^4~4b . 

Now we assess cluster lens-distor tion profiles from the 
color-selected background galaxies ( jj4.3b for the four clus- 
ters, in order to examine the form of the underlying cluster 
mass profile and to characterize cluster mass properties. In 
the weak lensing limit (n, \ j\ <C 1), the azimuthally averaged 
tangential distortion profile (g+(9)) (eq. 126*1) is related to the 
proje cted mass density profile (e.g., Bartel mann & Schneiderl 
2001) as 



(g + (6))~( 1+ (e)) = R«6)-(K(0)), 



(29) 



where (•••) denotes the azimuthal average, and k(< 9) is 
the mean convergence within a circular aperture of radius 
9 defined as R(< 9) = (tt^ 2 )" 1 lff, <e (P0' k(0'). Note that 
equation d29l holds for an arbitrary mass distribution. With 
the assumption of quasi-circular symmetry in the projected 
mass distribution, one can express the tangential distortion as 
(g+(9)) ~ [k(< 8)- (k(0))]/[1 - (k(0))] in the non-linear but 
sub-critical (det4(0) > 0) regime. 

Figure [5] shows the azimuthally-averaged radial profiles of 
the tangential distortion, (g+) (E mode), and the 45° -rotated 
component, (g x ) (B mode). Here the presence of B modes can 
be used to check for systematic errors. For each of the clus- 
ters, the observed E-mode signal is significant at the 12-16er 
level out to the limit of our data (9 ~ 20'). The significance 
level of the B-mode detection is about 2.5c for each cluster, 
which is about a factor of 5 smaller than E-mode. 
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The measured g+ profiles are compared with two represen- 
tative cluster mass models, namely the NFW model and the 
singular isothermal sphere (SIS) model. Firstly, the NFW uni- 
versal density profile has a two-parameter functional form as 



PNFWO) = 



( r /r s )(l + r/r s ) 2 



(30) 



where p s is a characteristic inner density, and r s is a character- 
istic inner radius. The logarithmic gradient n = dlnp(r)/dlnr 
of the NFW density profile flattens continuously towards the 
center of mass, with a flatter central slope n = — 1 and a steeper 
outer slope (n — ► -3 when r — > oo) than a purely isother- 
mal structure (n = -2). A useful index, the concentration, 
compares the virial radius, r v i r , to r s of the NFW profile, 
Cvk = fvir/rj. We specify the NFW model with the halo virial 
mass M v ir and the concentration c yu instead of p s and r s . is 
We employ the radial dependence of the NFW lensing pro- 
files, knfw(^) an d 7+ nfw( 0), given bv lBartermannld 19961) and 
IWright & Brainerdl (l2000h . Next, the SIS density profile is 
given by 



psi${r) = 



2nGr 2 



(31) 



where er v is the one-dimensional isothermal velocity disper- 
sion of the SIS halo. The lensing profiles for the SIS model, 
obtained by projections of the three-dimensional mass distri- 
bution, are found to be 



ksis(#) = 7+,sis(60 = — , 



(32) 



where 9e is the Einstein radius defined by 9e = 
4iT(a v /c) 2 D ds /D s . 

Table [3] lists the best-fitting parameters for these models, 
together with the predicted Einstein radius 9e for a fiducial 
source at z s = 1.5, corresponding roughly to the median red- 
shifts of our blue background galaxies. For a quantitative 
comparison of the models, we introduce as a measure of 
the goodness-of-fit the significance probability Q(u/2, x 2 /2) 
to find by chance a value of \ 2 as poor as th e observed 
value for a given number of dof, v (see §15.2 in iPress et al.l 
I1992I) . 19 We find with our best-fit NFW models Q- values of 
Q ~ 0.50,0.95,0.36, and 0.80, and with our best-fit SIS mod- 
els Q ~ 0.28,5.0 x 10" 6 ,0.37, and 0.87, for A1689, A2142, 
A2261, and A2390, respectively. Both models provide sta- 
tistically acceptable fits for A 1689, A2261, and A2390. For 
our lowest-z cluster A2142, the curvature in the g + profile is 
pronounced, and a SIS model for A2142 is strongly ruled out 
by the Subaru distortion data alone, where the minimum % 2 is 
XL = 39 with 8 dof. 

4.5.2. Lens Convergence 

Although the lensing distortion is directly observable, the 
effect of substructure on the gravitational shear is non- 
local. Here we examine the lens convergence (k) profiles us- 
ing the shear-based ID recon struction metho d dev eloped by 
Umetsu & Broadhurst (2008). See Appendix lA.ll for details 
of the reconstruction method. 

18 We assume the cluster redshift Zd is equal to the cluster virial redshift. 

19 Note that a Q value greater than 0.1 indicates a satisfactory agreement 
between the data and the model; if Q it, 0.001, the fit may be acceptable, 
e.g. in a case that the measurement errors are moderately underestimated; if 
Q ^ 0.001, the model may be called into question. 



In Figure [6] we show, for the four clusters, model- 
independent k profiles derived using the shear-based ID re- 
construction method, together with predictions from the best- 
fit NFW models for the k(6) and g+(9) data. The substructure 
contribution to k(9) is local, whereas the inversion from the 
observable distortion to k involves a non-local process. Con- 
sequently the ID inversion method requires a boundary con- 
dition specified in terms of the mean k value within an outer 
annular region (lying out to 18'-19')- We determine this value 
for each cluster using the best-fit NFW model for the g + pro- 
file (Table 0]>. 

We find that the two sets of best-fit NFW parameters are 
in excellent agreement for all except A226 1 : For A226 1 , the 
best-fit values of c v i r from the g+ and k profiles are are in 
poorer agreement. From Figures [4] and [6] we see that the 
NFW fit to the g+ profile of A2261 is affected by the pres- 
ence of mass structures at outer radii, 6 ~ 4' and 10', result- 
ing in a slightly shallower profile (c v ; r ~ 6.4) than in the k 
analysis. It turns out that these mass structures are associ- 
ated with galaxy overdensities whose mean colors are redder 
than the cluster sequence for A2261 at z = 0.224, A(V-R C ) = 
(V - R c ) - (V - R c )a226i ~ +0.6, and hence they are likely to 
be physically unassociated background objects. The NFW 
fit to k(9) yields a steeper profile with a high concentration, 
c v j r ~ 10.2, which implies a large Einstein angle of 9% ~ 37" 
at z s = 1 -5 (TableO. This is in good agreement with our pre- 
liminary strong-lensing model (Zitrin et al. in p reparation) 
based on the method by Broadhur st et al.l d2005al) . in which 
the deflection field is constructed based on the smoothed clus- 
ter light distribution to predict the appearance and positions 
of counter images of background galaxies. This model is re- 
fined as new multiply-lensed images are identified in deep 
Subaru VR C and CFHTAVIRCam JHK S images, and incorpo- 
rated to improve the cluster mass model. Figure [7] shows the 
tangential critical curve predicted for a background source at 
z s ~ 1 .5, overlaid on the Subaru V +R C pseudo-color image in 
the central 6.7' x 6.7' region of A2261. The predicted critical 
curve is a nearly circular Einstein ring, characterized by a n 
effective radius of 9 E ~ 40" (see lOguri & Blandfordll2008l) . 
This motivates us to further improve the statistical constraints 
on the NFW model by combining the outer lens convergence 
profile with the observed constraint on the inner Einstein ra- 
dius. A joint fit of the NFW profile to the n profile and the 
inner Einstein-radius constraint with 6% = 40" ±4" (z s =1) 
tightens the constraints on the NFW parameters (see §5.4.2 of 
Umetsu & Broadhurst 2008): M vir = 1.25^;^ x IO^Mq/T 1 
and c v ; r = ll.l*i 9> This model yields an Einstein radius of 
9e = (40 ± 1 1)" at z s = 1 .5. In the following analysis we will 
adopt this as our primary mass model of A2261. 

For the strong-lensing cluster A 1689, more detailed 
lensing constraints are available from joint observa- 
tions with the high-resolution Hubble Space Telescope 
(HST) Advanced Camera for Surveys (ACS) and the 
wide-field Subaru/Suprime-C am dBrpadhurst et alj fe OOSb; 
lUmetsu & Broadhurstl 120081) . In lUmetsu & Broadhursfl 
(2008) we combined all possible lensin g measurements, 
namely , the ACS strong-lensing profile of Broadhur st et alj 
(2005b) and Subaru weak-lensing distortion and magnifica- 
tion data, in a full two-dimensional treatment, to achieve the 
maximum possible lensing precision. Note, the combination 
of distortion and magnification data breaks the mass-sheet de- 
generacy (see eq. |[T6l ) inherent in all rec onstruction meth- 
ods based on distortion information alone (Bartelmann et al. 
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1996). It was found that the joint ACS and Subaru data, cov- 
ering a wide range of radii from 10 up to 2000kpc/T I , are 
well approximated by a single NFW profile with M v i r = (1 .5 ± 
O.l^f) x IO^Mq/T 1 and c vir = 12.7 ± 1 ±2.8 (statistical fol- 
lowed by systematic uncertainty at 68% confidence). 20 This 
properly reproduces the Einstein radius, whi ch is tightly con- 
strained by detailed strong-lens modeling dBro adhur st et all 
l2005at lHalkola et al.ll2006t iLimousin et al.ll2007l) : 6v ~ 52" 
at z s = 3.05 (or 8e — 45" at a fiducial source redshift of 
z s = 1). With the improved color selection for the red back- 



ground sample (see 134.31 ). we have redone a joint fit to the 
A CS and Subaru lensing obser vations using the 2D method 
of lUmetsu & Broadhurst (2008): The refined constraints on 
the NFW parameters are M vir = 1.55^;[| x 1O 15 M /i _1 and 
c v ; r = 12. 3^] g, yielding an Einstein radius of 50^q arcsec at 
Z s = 1-5. In the following, we will adopt this refined NFW 
profile as our primary mass model of A1689. 

5. DISTRIBUTIONS OF MASS AND HOT BARYONS 

Here we aim to compare the projected distribution of mass 
and ICM in the clusters using our Subaru weak lensing and 
AMiBA SZE maps. To make a quantitative comparison, we 
first define the "cluster shapes" on weak lensing mass struc- 
ture by introducing a spin-2 halo ellipticity e hal ° = + z'eij °> 
defined in terms of weighted quadrupole shape moments 
(a, (3= 1,2), as 



„halo 



(M = 



halo 



Qn 



ehalo 
22 



halo • 



ehalo I /Ohal 
11 + ^22 



ehalo I /ohalo 
1 1 + «22 



(33) 



Q^°(^p)= / d 2 9A9 a A6f3 k(0), (34) 
JAe<e il? 

where 6* ap is the circular aperture radius, and A8 a is the angu- 
lar displacement vector from the cluster center. Similarly, the 
spin-2 halo ellipticity for the SZE is defined using the cleaned 
SZE decrement map -Al{6) cx y(9) instead of n{0) in equa- 
tion (1341 . The degree of halo ellipticity is quantified by the 

modulus of the spin-2 ellipticity, |e halo | = y / (e halo ) 2 + (e5 al °) 2 , 
and the orientation of halo is defined by the position angle 
of the major axis, hal ° = arctan(e2 alo /<? 1 j lalo )/2. In order to 
avoid noisy shape measurements, we introduce a lower limit 
of k(8) > and -Al(0) > in equation d34l . Practical shape 
measurements are done using pixelized lensing and SZE maps 
shown in Figure |4] The images are sufficiently oversampled 
that the integral in equation d34i l can be approximated by the 
discrete sum. Note, a comparison in terms of the shape param- 
eters is optimal for the present case where the paired AMiBA 
and weak lensing images have different angular resolutions: 
6Whm - 6' FWHM for AMiBA7, and 6Whm ^ 1 -7' FWHM 
for Subaru weak lensing. When the aperture diameter is larger 
than the resolution #fwhm, i.e., #ap > #fwhm/2, the halo shape 
parameters can be reasonably defined and measured from the 
maps. 

Now we measure as a function of aperture radius # ap 
the cluster ellipticity and orientation profiles for projected 
mass and ICM pressure as represented by the lensing k 
and SZE decrement maps, respectively. For the Subaru 
weak lensing, the shape parameters are measured at ap = 



20 In Umetsu & Broadhurst 12008) cluster masses are expressed in units 
of 10 15 Mq with h = 0.7. The systematic uncertainty in M v j, is tightly cor- 
related with that in c v j r through the Einstein radius constraint by the ACS 
observations. 



[1,2,3,..., 11] x 6Whm (l-7'£0ap £ 18.3'); for the AMiBA 
SZE, ap = [1,2,3] X 6Whm (6' £ 6> ap £ 18'). The level of un- 
certainty in the halo shape parameters is assessed by a Monte- 
Carlo error analysis assuming Gaussian errors for weak lens- 
ing distortion and AMiBA visibili ty measurements (f or the 
Gaussianity of AMiBA data, see iNishioka et al. 1 120081) . For 
each cluster and dataset, we generate a set of 500 Monte Carlo 
simulations of Gaussian noise and propagate into final uncer- 
tainties in the spin-2 halo ellipticity, e hal °. Figure [8] displays, 
for the four clusters, the resulting cluster ellipticity and orien- 
tation profiles on mass and ICM structure as measured from 
the Subaru weak lensing and AMiBA SZE maps, shown sep- 
arately for the ellipticity modulus |<? halo | and the orientation, 
2^,haio (]- w j ce the position angle). Overall, a good agreement 
is found between the shapes of mass and ICM structure up to 
large radii, in terms of both ellipticity and orientation. In par- 
ticular, our results on A2142 and A2390 show that the mass 
and pressure distributions trace each other well at all radii. 
At a large radius of £? ap ^ 10', the position angle of A2142 is 
^haio _ 5Q o For A2 390, the position angle is hal ° - 30° at 
all radii. 

6. CLUSTER GAS MASS FRACTION PROFILES 

6.1. Method 

In modeling the clusters, we consider two representative 
analytic models for descr ibing the cluster DM and ICM dis- 
tributions, namely (1) the lKomatsu & Seliakl d2001l hereafter 
KS01) model of the universal gas density /temperature profiles 
and (2) the isothermal (3 model, where both are physically mo- 
tivated under the hypothesis of hydrostatic equilibrium and 
polytropic equation-of-state, P cx p 1 , with an additional as- 
sumption about the spherical symmetry of the system. 

Joint AMiBA SZE and Subaru weak lensing observations 
probe cluster structures on angular scales up to AO ~ 23'. 21 
At the median redshift z — 0.2 of our clusters, this maxi- 
mum angle covered by the data corresponds roughly to r2oo ~ 
0.8r v i r , except rsoo ~ 0.5r v ; r for A2142 at z = 0.09. In order 
to better constrain the gas mass fraction in the outer parts 
of the clusters, we adopt a prior that the gas density pro- 
file p gas (r) traces the underlying (total) mass density profile, 
PtotM- Such a relationship is expected at large radii, where 
non-gravitational processes, such as radiative cooling and star 
formation, have not had a major effect on the structure of the 
atmosphere so that the polytropic assumption remains valid 
(Lewis et al. 2000). Clearly this results in the gas mass frac- 
tion, p g as(r)/ Ptoi(r), tending to a constant at large radius. In 
the context of the isothermal (3 model, this simply means that 
= 2/3. 

In both models, for each cluster, the mass density profile 
PtotM is constrained solely by the Subaru weak lensing data 
(®, the gas temperature profile r gas (r) is normalized by the 
spatially-averaged X-ray temperature (see Table [TJ, and the 
electron pressure profile P e (r) = n e (r)lcBT e (r) is normalized by 
the AMiBA SZE data, where n e (r) is the electron number den- 
sity, and T e (r) = T gas (r) is the electron temperature. The gas 
density is then given by p gas (r) = p, e m p n e (r). 

6.2. Cluster Models 
6.2.1. NFW-Consistent Model of Komatsu & Seljak 2001 

21 The FWHM of the primary beam patten of the AMiBA is about 23', 
while the field-of-view of the Subaru/Suprime-Cam is about 34' 
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The KS01 model describes the poly tropic gas in hydro- 
static equilibrium with a gravitational potential described 
by the universal density profile for c ollisionless CDM 
halos prop osed by iNavarro et al] (| 19961 hereafter NFW). 
See KS01, [Komatsu & Seliakl fl2002l hereafter KS02), and 
IWorrall &~B irkinshaw (2006) for more detailed discussions. 
High mass clusters with virial masses M v j r £ 10 15 M Q //z are 
so massive that the virial temperature of the gas is too high 
for efficient cooling and hence the cluster potential simply re- 
flects the dominant DM. This has been recently established 
by our Subaru weak lensing study of several massive cluster s 
dBroadhurstet al.ll2005Ul2008tlUmetsu & Broadhurstl2 008). 

In this model, the gas mass profile traces the NFW pro- 
file in the outer region of the halo (r vu /2 £ r ^ r vlr ; see 
KS01), satisfying the adopted prior of the constant gas 
mass fraction p sas (r)/ p tot (r) at large radii. This behav- 
ior is supp orted by cosmologica l hydrodynamic simula- 
tions (e.g., lYoshikawa et al.l 120001) . and is recently found 
from the stacked SZE analy sis of the WMAP 3 -year data 
(lAtrio-Barandela et al.ll2008h . The shape of the gas distribu- 
tion functions, as well as the poly tropic index 7 gas , can be 
fully specified by the halo virial mass, M v ; r , and the halo con- 
centration, c v ; r = r v ; r / r s , of the NFW profile. 

In the following, we use the form of the NFW profile to 
determine r wu , r2oo, '"SOO; and r2soo- Table [4] summarizes the 
NFW model paramete rs de rived from our lensing analysis for 
the four clusters (see 34.51 l. For each cluster we also list the 



corresponding (r 25 oo,?"500,?"200,?"vir)- For calculating 7 gas and 
the normalization factor 7/(0) for a structure constant (B in 
equation [16] of KS02), we follow the fitting formulae given 
by KS02, which are valid for halo concentration, 1 < c vir < 25 
(see Table |4j. For our clusters, 7 gas is in the range of 1.15 
to 1.20. Following the prescription in KS01, we convert the 
X-ray cluster temperature Tx to the central gas temperature 
7g as (0) of the KS01 model. 

6.2.2. Isothermal (3 Profile 

The isothermal (3 model provides an alternati ve consistent 
soluti on of the hydrostatic equilibrium equation (Hatt oriet al.l 
fT999h . assuming the ICM is isothermal and its density pro- 
file follows p gas (r) = /3 gas (0)[l + (r/r c ) 2 ]~ 3 P/ 2 with the gas core 
radius r c . At large radii, r ^> r c , where both of our SZE and 
weak lensing observations are sensitive, the total mass density 
follows ptot(r) oc r~ 2 . Thus we set (3 = 2/3 to satisfy our as- 
sumption of constant p g - ds (r) / p tot (r) at large radius. We adopt 
the values of r c and Tx listed in Table Q] taken from X-ray 
observations, and use r gas (r) = Tx as the gas temperature for 
this model. At r 3> r c , the pt o t (r) profile can be approximated 
by that of a SIS (see 34.5. U parametrized by the isothermal 
ID velocity dispersion <j v (see Table 0J, constrained by the 
Subaru distortion data (see 34.5b . 

Requiring hydrostatic balance gives an isothermal temper- 
ature Tsis, equivalent to <r v , as 

2 2 

k B T S i S = pm p a v — . (35) 

For j3 = 2/3, fcgTsis = pm p a 2 ,, which can be compared with 
the observed Tx (Table [TJ. For our AMiBA-lensing cluster 
sample, we found X-ray to SIS temperature ratios Tx/Tsis = 
0.82 ± 0.03, 1 .65 ± 0.15,0.94 ± 0.05, 1 .28 ± 0.15 for A1689, 
A2142, A2261, and A2390, respectively. For A2261 and 
A2390, the inferred temperature ratios are consistent with 
unity at I— 2a. For the merging cluster A2142, the ob- 
served spatially-averaged X-ray temperature (cooling-flow 



corrected; see Mark evitch|[T998l) is significantly higher than 
the lensing-derived temperature. This temperature excess of 
~ 4er could be explained by the eff ects of merger boosts, as 
discussed in Okabe & Umetsul d2008l) . The temperature ratio 
Tx/Tsis for A1689, on the other hand, is significantly lower 
than uni ty. Recently, a simila r level of discrepancy was also 
found in lLemze et al.l d2008al) . who performed a careful joint 
X-ray and lensing analysis of this cluster. A deprojected 3D 
temperature profile was obtained using a model-independent 
approach to the Chandra X-ray emission measurements and 
the projected mass profile obtained from the j oint strong/weak 
lensing analysis of Broadhurs t et al.l (l2005bl) . The projected 
temperature profile predicted from their joint analysis ex- 
ceeds the observed temperature by 30% at all radii, a level of 
discrepancy suggested from hydrodynamical simulations that 
find that denser, more X-ray luminous small-scale structure 
can bias X-ray t emperature measure ments downward at about 
the same level (l^wahara et al. 2007]). If we accept this +30% 
correction for T x , the ratio T X /T SIS -> 1 .07 ± 0.04 for A1689, 
consistent with (3 = 2/3. 

6.3. AM IB A SZE Data 

We use our AMiBA data to constrain the remaining nor- 
malization parameter for the p g - ds {r) profile, p gas (0). The cal- 
ibrated output of the AMiB A interferomete r, after the lag- 
to-visibility transformation dWu et a l. 2008a), is the complex 
visibility V(u) as a function of baseline vector in units of 
wavelength, u = d/X, given as the Fourier transform of the 
sky brightness distribution Al(6) attenuated by the antenna 
primary beam pattern A(9). 

In targeted AMiBA observations at 94 GHz, the sky signal 
Al(6) with respect to the background (i.e., atmosphere and 
the mean CMB intensity) is dominated by the thermal SZE 
due to hot electrons in the cluster, A/sze = ^norm^(^)y (see 
eq. H]). The Comptonization parameter y is expressed as 
a line-of-sight integral of the thermal electron pressure (see 
eq. J2]). In the line-of-sight projection of equation (O, the 
cutoff radius r max needs to be specified. We take r max = a r r v ; r 
with a dimensionless constant a r which we set to a r = 2. 
In the present study we found the line-of-sight projection 
in equation (f2]l is insensitive to the choice of a r as long as 
a r >l. 

A useful measure of the thermal SZE is the integrated 
Comptonization parameter Y(8), 

Y(6) = 2tt [ de'e'y(e'), (36) 
Jo 

which is proportional to the SZE flux, and is a measure of 
the thermal energy content in the ICM. The value of Y is 
less sensitive to the details of the model fitted than the central 
Comptonization parameter yo = y(0), with the current config- 
uration of AMiBA. If the A(0)y(6) field has reflection sym- 
metry about the pointing center, then the imaginary part of 
V(u) vanishes, and the sky signal is entirely contained in the 
real visibility flux. If the A(9)y(9) field is further azimuthally 
symmetric, the real visibility flux is expressed by the Hankel 
transform of order zero as 

/Re/ 



V Ke (u) = 2irl nmm g(v c ) / d9 6A{6)y(6)J () (2Ku6) 



i2irl 



d98A(6)—M2iru9), 

yo 



(37) 



where Iq = Inoimg(y c )yo is the central SZE intensity at v c = 
94 GHz, Jq(x) is the Bessel function of the first kind and or- 
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der zero, and A(8) is well approximated by a circular Gaus- 
sian with FWHM = 1.22(A/D) ^ 23' at v c = 94G Hz with an 
antenna diameter of Z) = 60cm dWu et alj|2008al) . The ob- 
served imaginary flux can be used to check for t he effects 
of pri mary CMB and radio source contamination dLiu et al.1 
2008). From our AMiBA data we derive in the Fourier do- 
main azimuthally-averaged visibility profiles (V(u)) for indi- 
vidual clusters. 

We constrain the normaliza tion Iq from \ 2 fitting to the 
(V(u)) profile dLiu et al.ll2008l) . In order to convert Iq into 
the central Comptonization parameter, we take account of 
(i) the relativistic correction 5sze(v, r gas ) in the SZE spec- 
tral function g(v) (see eq. Q) and (ii) cor rections for con - 
tamination by discrete radio point sources dLiuetalJl2008h . 
The level of contamination in Iq from known discrete point 
sources has been estimat ed to be about (6-35)% in our four 
clusters dLiu et al.ll20 08). In all cases, a net positive contri- 
bution of point sources was found in our 2-patch differencing 
AMiBA observations (§0, indicating that there are more ra- 
dio so urces towards clusters than in the background dLiu et al.l 
2008). Thus ignoring the point source correction would sys- 
tematically bias the SZE flux estimates, leading to an under- 
estimate of yo. The relativistic correction to the thermal SZE 
is 6-7% in our Tx range at 94 GHz. 

Table [5] summarizes, for our two models, the best-fitting 
parameter, yo, and the T -parameter interior to a cylinder of 
radius 9 = 3' that roughly matches the AMiBA synthesized 
beam with 6' FWHM. For each case, both cluster models yield 
consistent values of yo and Y(3') within ler; in particular, the 
inferred values of Y(3') for the two mod els a re in excellent 
agreement. Following the procedure in $6.11 we convert yo 
into the central gas mass density, /O gas (0). 

6.4. Gas Mass Fraction Profiles 



We derive cumulative gas fraction profiles, 



/gas(< r) ■■ 



M gas (< r) 
M to t« r) ' 



(38) 



for our cluster sample using two sets of cluster models de- 
scribed in $6.21 where M gas (< r) and M tot (< r) are the hot gas 
and total cluster masses contained within a spherical radius r. 
In Table|6]we list, for each of the clusters, M gas and / gas within 
7"2500, r 500, and r^ao (see also Table|4]l calculated with the two 
models. Note that our total mass estimates do not require the 
assumption of hydrostatic balance, but are determined based 
solely on the weak lensing measurements. Gaussian error 
propagation was used to derive the errors on M gas (< r) and 
/gas(< r )- We propagate errors on the individual cluster pa- 
rameters (Tables Q] and |4]i by a Monte-Carlo method. For 
A2142, the isothermal model increasingly overpredicts / gas 
at all radii r > r25oo, exceed ing the cosmic bary on fraction 
ft = n b /Q, m = 0.171 ±0.009 dDunklev et alJ fcOOS). For other 
clusters in our sample, both cluster models yield consistent 
/ gas and M gas measurements within the statistical uncertain- 
ties from the SZE and weak lensing data. 

Our SZE/weak lensing-based measurements can be com- 
pared w ith other X-ray and SZE measurements. Gre go et ail 
(2001b) derived gas fractions for a sample of 18 clus- 
ters from 30 GHz SZE observations with BIMA and OVRO 
in combination with published X-ray temperatures. They 
found / gas « r 500 ) = 0.140t<*% and 0.053^/^ (h 10 = 
h/0.7) for A1689 and A2261, respectively, in agreement with 
our results. For A2142, the / gas and M gas values inferred 
from the KS01 model are in good agreement with those 



from the VSA SZE observations at 30 GHz (Lancas ter et alj 
HH, M gas (r 50 o) = 6.l + _\i x 1O 13 M /T 2 and / gas (r 5 oo) = 
0.12 3^n^7n- From a d etailed analysis of Chandra X-ray 
data, IVikhlinin etaT] d2006l) obtained / gas (< r 500 ) = (0.141 ± 

0.009)/i72 /2 (h 72 = h /0.72) for A2390, in good agreement with 
our results. 

Furthermore, it is interesting to compare our results with 
the detailed j oint X- ray and lensing analysis of A1689 by 
iLemze et al.1 d2008al) . who derived deprojected profiles of 
Aot(r), /9 gas (r), and T gRS (r) assuming hydrostatic equilibrium, 
using a model-independent approach to the Chandra X-ray 
emission profile and the projected lensing mass profile of 
iBroadh urst et al. (2005b). A steep 3D mass profile was ob- 
tained by this approach, with the inferred concentration of 
c v jr = 12.2^? n, consistent with the detailed lensing analysis o f 
IBroadhurst etai] d2005bl) and lUmetsu & Broadhurstl d2008l) . 
whereas the observed X-ray temperature profile falls short of 
the derived profile at all radii by a constant factor of ~ 30% 
(see 36.2.2| i. With the pressure profile of Lemze et al. we 
find yo = (4.7 ± 0.3) x 10~ 4 , which is in agreement with our 
KS01 prediction, y = (4.2 ± 1.0) x 10" 4 (Table 0. The inte- 
grated Comptonization parameter predicted by the Lemze et 
al. modelis Y(3') = (3.0±0.1) x 10" 10 , which roughly agrees 
with the AMiBA measurement of 7(3') = (2.5 ± 0.6) x 10" 10 . 
Alternatively, adopting the observed temperature profile in 
the Lemze et al. model reduces the predicted SZE signal 
by a factor of - 30%, yielding y ~ 3.3 x 10" 4 and Y(3') ~ 
2.1 x 10~ 10 , again in agreement with the AMiBA measure- 
ments. Therefore, more accurate SZE measurements are re- 
quired to further test and verify this detailed cluster model. 

We now use our data to find the average gas fraction pro- 
file over our sample of four hot X-ray clusters. The weighted 
average of M v ; r in our AMiBA-lensing sample is (M v i r ) = 
(1.19 ±0.08) x l0 15 M Q h~ l , with a weighted-mean concen- 
tration of (c V ir) = 8.9 ±0.6. The weighted average of the 
cluster virial radius is (r v - n ) ~ 1.95MpcA -1 . At each radius 
we compute the sample-averaged gas fraction, (/ gas (< r)), 
weighted by the inverse square of the statistical ler uncer- 
tainty. In Figure [9] we show for the two models the resulting 
(/gas) profiles as a function of radius in units of r v ; r , along 
with the published results for other X-ray and SZE obser- 
vations. Here the uncertainties (cross-hatched) represent the 
standard error (ler) of the weighted mean at each radius point, 
including both the statistical measurement uncertainties and 
cluster-to-cluster variance. Note A2 142 has been excluded for 
the isothermal case (see above). The averaged (/ gas ) profiles 
derived for the isothermal and KS01 models are consistent 
within lc at all radii, and lie below the cosmic baryon frac- 
tion/, = 0.171 ±0.00 9 constrained by the WMAP 5-year data 
dDunklev etalJl2008h . At r = (r 20 o) ^ 0.79(r vir ), the KS01 
model gives 



(/ g as,20o) = 0.133 ±0.020±0.018 



(39) 



where the first error is statistical, and the second is the 
standard error due to cluster-to-cluster variance. This is 
marginally consistent with (/ gas ,20o) = 0.109 ±0.013 ob- 
tained from the averaged SZE profile of a sample of 193 
X-ray clusters (Tx > 3keV) using the WMAP 3 -year data 
dAfshordi et al.ll2007l) . A similar value of (/ gas .2 oo) = 0.11 ± 
0.03 was obtained by iBiviano & Saluccil d2006l) for a sam- 
ple of 59 nearby clusters from the ESO Nearby Abell Clus- 
ter Survey, where the total and ICM mass profile are deter- 
mined by their dynamical and X-ray analyses, respectively. 
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At r = (rsoo) ~ 0.53 (r v i r ), we have 

</ gas ,50o) = 0. 126 ± 0.019 ± 0.016 (40) 

for the KS01 model, in good agreement with the Chandra 
X -ray measurements for a subset of six Tx > 5keV clusters 
in IVikhlinin et al l (120061) . At r = (r 25 oo) ^ 0.25 (r vir ), which 
is close to the resolution limit of AMiBA7, we have for the 
KS01 model 

(/gas,250o) =0.105 ±0.015 ±0.012, (41) 

again marginally consistent with the Chandra gas fraction 
measurements in 2 6 X-ray luminous clusters with Tx > 5 keV 
(lAllen et alJl2004b . 

7. DISCUSSION AND CONCLUSIONS 

We have obtained secure, model-independent profiles of 
the lens distortion and projected mass (Figures [5] and |6j by 
using the shape distortion measurements from high-quality 
Subaru imaging, for our AMiBA lensing sample of four high- 
mass clusters. We utilized weak lensing dilution in deep Sub- 
aru color images to define color-magnitude boundaries for 
blue/red galaxy samples, where a reliable weak lensing sig- 
nal can be derived, free of unlensed cluster members (Figure 
0. Cluster contamination otherwise preferentially dilutes the 
inner lensing signal leading to spuriously shallower profiles. 
With the observed lensing profiles we have examined cluster 
mass-density profiles dominated by DM. For all of the clus- 
ters in our sample, the lensing profiles are well described by 
the NFW profile predicted for collisionless CDM halos. 

A qualitative comparison between our weak lensing and 
SZE data, on scales r £ 3' limited by the current AMiBA res- 
olution, shows a good correlation between the distribution of 
mass (weak lensing) and hot baryons (SZE) in massive clus- 
ter environments ( §I4.41 i. as physically expected for high mass 
clusters with deep gravitational potential wells Q4.41 >. We 
have also examined and compared, for the first time, the clus- 
ter ellipticity and orientation profiles on mass and ICM struc- 
ture in the Subaru weak lensing and AMiBA SZE observa- 
tions, respectively. For all of the four clusters, the mass and 
ICM distribution shapes are in good agreement at all relevant 
radii in terms of both ellipticity and orientation (Figure|H). In 
the context of the CDM model, the mass density, dominated 
by collisionless DM, is expected to have a more irregular and 
elliptical distribution than the ICM density due to inherent tri- 
axiality of CDM halos. We do not see such a tendency in our 
lensing and SZE datasets, although our ability to find such ef- 
fects is limited by the resolution of the current AMiBA SZE 
measurements. 

We have obtained cluster gas fraction profiles (Figure 
[9]l for the AMiBA-lensing sample (T x > 8keV) based on 
joint AMiBA SZE and Subaru weak lensing observations 
( jj6.41 >. Our cluster gas fraction measurements are overall 
in good agreement with previously-published values. At 
r = (f"20o) — 0.79 (r v ir), corresponding roughly to the maxi- 
mum available radius in our joint SZE/weak lensing data, 
the sample-averaged gas fraction is (/ gas ,20o) = 0.133 ±0.027 
for the NFW-consistent KS01 model, representing the aver- 
age over our high-mass cluster sample with a mean virial 
mass of (M vir ) = (1.2±0.1) x 1O 15 M A -1 . When compared 
to the cosmic baryon fraction ft,, this indicates (/ gas ,20o) //* = 
0.78 ± 0.16, i.e., (22 ± 16)% of the bar yons is missing from 
the hot phase in our cluster sample (cf. Afshordi et al.ll2007t 
ICrain et al.ll2007h . This missing cluster baryon fraction is par- 
tially made up by observe d stellar and cold gas f ractions of ~ 
several % in our Tx range (Gonzalez et al. 2005). 



Halo triaxiality may affect our projected total and gas mass 
measurements based on the assumption of spherical symme- 
try, producing an orientation bias. A degree of triaxiality is in- 
evitable for collisionless gravitationally collapsed structures. 
The likely effect of triaxiality on the measure ments of lensing 
properties has been examined analy tically dQguri et al1 l2005b 
ISerenoll2007t ICorless & K ing 2007 ]), and in numer ical inves- 
tigations dJing & Sutol l2002t iHennawi et al.ll2007l) . The ef- 
fect of triaxiality will be less for the collisional ICM, which 
follows the gravitational potential and will be more spheri- 
cal and more smoothly distributed than the total mass density 
distribution. For an unbiased measurement of the gas mass 
fractions, a large, homogeneous sample of clusters would be 
needed to beat down the orientation bias. 

Possible biases in X-ray spectroscopic temperature mea- 
surements (Maz zotta et al.ll2004t iKawahara et al.ll2007h may 
also affect our gas fraction measurements based on the over- 
all normalization by the observed X-ray temperature. This 
would need to be taken seriously into account in future inves- 
tigations with larger samples and higher statistical precision. 

Our joint analysis of high quality Subaru weak lensing 
and AMiBA SZE observations allows for a detailed study 
of individual cluste rs. The cluster A2142 shows complex 
mass substructure (Okab e~& Umetsul 12008). and displays a 
shallower density profile with c wu ~ 5, consistent with de- 
tailed X-ray observations which imply recent interaction. Due 
to its low-z and low c vlI , the curvature in the lensing pro- 
files is highly pronounced, so that a SIS profile for A2142 
is strongly ruled out by the Subaru distortion data alone 
( 34.5. 11 1. For this cluster, our AMiBA SZE map shows an 
extended structure in the ICM distribution, elongated along 
the northwest-southeast direction. This direction of elonga- 
tion in the SZE halo is in good agree ment with the cometary 
X-ray appearance seen by Chandra dMarkevitch et al]l2000t 
lOkabe & Umetsul 120081) . In addition, an extended structure 
showing some excess SZE can be seen in the northwest region 
of the cluster. A joint weak-lensing, opti cal-photometric, and 
X-ray analysis ( lOkabe & Umetsul 120081) revealed northwest 
mass substructure in this SZE excess region, located ahead 
of the northwest edge of the central gas core seen in X-rays. 
The northwest mass substructure is also seen in our weak lens- 
ing mass map (Figure |4]i based on the much improved color 
selection for the background sample. A slight excess of clus- 
ter sequence galaxie s associated with t he nor thwest substruc- 
ture is also found in Oka be & Umetsul ((2008), while no X-ray 
counterpart is seen in the Chandra and XMM-Newton images 
(Okabe & Umetsu 2008). Good consistency found between 
the SZE and weak lensing maps is encouraging, and may sug- 
gest that the northwest excess SZE is a pressure increase in 
the ICM associated with the moving northwest substructure. 
Clearly further improvements in both sensitivity and resolu- 
tion are needed if SZE data are to attain a significant detec- 
tion of the excess structure in the northwest region. Nonethe- 
less, this demonstrates the potential of SZE observations as a 
powerful tool for measuring the distribution of ICM in cluster 
outskirts where the X-ray emission measure (ex rfy is rapidly 
decreasing. This also demonstrates the potential of AMiBA, 
and the power of multiwavelength cluster analysis for prob- 
ing the distribution of mass and baryons in clusters. A further 
detailed multiwavelength analysis of A2142 will be of great 
importance for further understanding of the cluster merger dy- 
namics and associated physical processes of the intracluster 
gas. 

For A2390 we obtain a highly elliptical mass distri- 
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bution at all radii from both weak and strong lensing 
dFrye & Broadhurstlll998l) . The elliptical mass distribution 
agrees well with the shape seen by AMiBA in the thermal 
SZE (Figures [4] and |8). Our joint lensing, SZE, and X-ray 
modeling leads to a relatively high gas mass fraction for this 
cluster, /gas. 500 ~ 0.15 for the NFW-consisten t case, which is 
in good agreement with careful X-ray work bv lVikhlinin et alj 

d2006h . / gas ,5oo = (0.141 ±0.009)/^2 /2 

We have refined for A 1689 the statistical constrai nts on 
the NFW mass model of lUmetsu & Broadhurstl (120081) . with 
our improved color selection for the red background sam- 
ple, where all possible lensing measurements are com- 
bined to achieve the maximum possible lensing precision, 
M vil - = (1.55lgj|) x IO^Mq/T 1 and c vir = 12.3^| (quoted 
are statistical errors at 68% confidence level), confirming 
again the high concentration found by detailed lensing work 
(Broadhurs t et alj|2005bt [ Halkola et al. 2006; ILimousin etall 
I2007t IUmetsu & Broadhurst 2008). The AMiBA SZE mea- 
surements at 94 GHz support the compact structure in the 
ICM distribution for this cluster (Figure |4). Recently, good 
consistency was found between high- quality multiwavelengt h 
datasets ava i lable f or this cluster dLemze et al.1 [2008a b). 
iLemze et al.1 (|2008a) performed a joint analysis of Chandra 
X-ray, ACS strong lensing, and Subaru weak lensing mea- 
surements, and derived an improved mass profile in a model- 
independent way. Their NFW fit to the derived mass profile 
yields a virial mass of M vil = (1 .58 ± 0.15) x 10 15 M Q /T 1 and 
a high concentration of c vrr = 12.2^' y, both of which are in 
excellent agreement with our full lensing constraints. More 
recently, Lemz e~etaT1 d2008b) further extended their multi- 
wavelength analysis to combine their X-ray/lensing measure- 
ments with two dynamical datasets from VLT/VTRMOS spec- 
troscopy and Subaru/Suprime-Cam imaging. Their joint lens- 
ing, X-ray, and dynamical analysis provides a tight constraint 
on the cluster virial mass: 1.5 < M Yil /(10 15 M Q h~ l ) < 1.65. 
Their purely dynamical analysis constrains the concentration 
parameter to be c vrr > 9.8 for A1689, in agreement with our 
indepe ndent lensing a nalysis and the joint X-ray/lensing anal- 
ysis of ILemze etall (l2008al) . We remark that NFW fits to 
the Subaru outer profiles alone give consistent but somewhat 
higher concentrations, c v i r ~ 15 (Table [3j see also Umetsu & 
Broadhurst 2008 and Broadhurst et al. 2008). This slight 
discrepancy can be explained by the mass density slope at 
large radii (9 ^ 5') for A1689 being slightly steeper than the 
NFW pro file where the asymptotic d ecline tends to pnfw(^) oc 
r~ 3 (see iBroadhurst etall l2005bt iMedezinski et all 120071: 
Lemze et alj2008aF Umetsu & Broadhurstl2008l:lLemze et alj 
2008b). Recent detailed modeling by ISaxton & Wu l d2008l) 
suggests such a steeper outer density profile in stationary, self- 
gravitating halos composed of adiabatic DM and radiative gas 
components. For accurate measurements of the outermost 
lensing profile, a wider optical/near-infrared wavelength cov- 
erage is required to improve the contamination-free selection 
of background galaxies, including blue background galaxies, 
behind this rich cluster. 

Our Subaru observations have established that A2261 is 
very similar to A 1689 in terms of both weak and strong 
lensing properties: Our preliminary strong lens modeling re- 
veals many tangential arcs and multiply-lensed images around 
A2261, with an effective Einstein radius 6% ~ 40" at z ~ 
1.5 (Figure [7]), which, when combined with our weak lens- 
ing measurements, implies a mass profile well fitted by an 
NFW model with a concentration c v ; r ~ 10, similar to A 1689 



( lUmetsu & Broa dhurst 2008), and considerably higher than 
theoretically expected for the standard ACDM model, where 
c v i r ~ 5 is predicted for the most massive relaxed clusters 
with M vir £ 10 15 M W dBullock et alj 1200 U iNeto et alj 120071: 
iDuffy et alJl2008h . 

Such a high concentration is also seen in several other 
well-studied massive clusters from care f ul lensing work 
dGavazzi et alj 120031: iKneibetalJ 120031: rB roadhur st et all 
| 2005bt ILimousin et al.l 120071: ILemze et ail [2008a; 
Broadhurst et al.l 120081) . The orientation bias due to 
halo triaxiality can potentially affect the projected lensing 
measurements, and hence the lensin g-based concentration 
measurement (e.g., Ogu ri et ai1l2005l) . A statistical bias in 
favor of prolate structure pointed to the observer is unavoid- 
able at some level, as this orientation boosts the projected 
surface mass density and hence the lensing signal. In the 
context of the ACDM model, this leads to an increase of 
~ 18% in the mean v alue of the lensing-based concentrations 
dHennawi etalJl2007h . A larger bias of ~ 30 up to 50% is 
expected for CDM halos selected by the presence of large 
gravit ational arcs dHennawi et al.l 120071 : lOguri & Blandfordl 
2008). Our cluster sample is identified by their being X- 
ray/SZE strong, with the added requirement of the availability 
of high-quality multi-band Subaru/Suprime-Cam imaging 
(see jj3.21 >. Hence, it is unlikely that the four clusters are all 
particularly triaxial with long axes pointing to the observer. 
Indeed, in the context of ACDM, the highly elliptical mass 
distribution of A2390 would suggest that its major axis is 
not far from the sky plane, and that its true concentration is 
higher than the projected measurement c v ; r ~ 7. 

A chance projection of structure along the line-of-sight may 
also influence the lensing-based cluster parameter determina- 
tion. It can locally boost the surface mass density, and hence 
can affect in a non-local manner (see eq. l29l ) the tangen- 
tial distortion measurement that is sensitive to the total in- 
terior mass, if this physically unassociated mass structure is 
contained within the measurement radius. For the determina- 
tion of the NFW concentration parameter, it can lead to either 
an under or over-estimate of the concentration depending on 
the apparent position of the projected structure with respect 
to the cluster center. When the projected structure is well 
isolated from the cluster center, one way to overcome this 
problem is to utilize the convergence profile to examine the 
cluster mass profile, by locally masking out the contri bution 
of known foreground/background structure (see 34.5.2l for the 
case of A2261). 

The ongoing upgrade of AMiBA to 13-elements with 1.2m 
antennas will improve its spatial resolution and dynamic 
range, and the 13-element AMiBA (AMiBA13) will be sen- 
sitive to structures on scales down to 2', matching the angu- 
lar scales probed by ground-based weak lensing observations 
dUmetsu et alj|2004h . For our initial target clusters, joint con- 
straints with AMiBA7 and AMiBA 13 data will complement 
the baseline coverage, which will further improve our mul- 
tiwavelength analysis of the relation between mass and hot 
baryons in the clusters. A joint analysis of complementary 
high-resolution lensing, SZE, and X-ray observations will be 
of great interest to address the issue of halo triaxiality and 
further impro ve the constraints on cluster density profiles 
dSerenoll2007l) . The AMiBA upgrade will also make the in- 
strument faster by a factor of ~ 60 in pointed observations. 
Our constraints can be further improved in the near future by 
observing a larger sample with AMiBA13. A detailed com- 
parison between X-ray based and SZE/weak lensing-based 
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gas fraction measurements will enable us to test the degree 
of dumpiness ((n^)/(n e ) 2 ) and of hydrostatic balance in hot 
cluster gas. The high angular resolution (2') of AMiBA13 
combined with dynamically-improved imaging capabilities 
will allow for direct tests o f the gas pressure p rofile in deep 
single pointed observations (Moln aret al.ll2008l) . 
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TABLE 1 

Target clusters and AMiBA/X-ray properties 



Cluster z 1 arcmin 3 
(kpc/T 1 ) 



A1689 0.183 129.6 

A2142 0.091 71.4 

A2261 0.224 151.8 

A2390 0.228 153.3 



AMiBA7 b 

SZE flux Image FWHM 

(mjy) (arcmin) 

-168 ±28 5.7 

-316 ±23 9.0 

-90 ±17 5.8 

-158 ±24 8.0 



X-ray c 

T x C Refs 

(keV) (arcmin) 



9.66^|o 0.44 ±0.01 3 

9.7±1.0 3.14±0.22 1,4,5 

8.82 + ^| 0.26 ±0.02 3 

10.1 ±1.1 0.47±0.05 2 



REF ERE NCES. — [1 ] Markevitc hetai] (19911); [2] Boehring eTeTaT] fi998l) : [3] IReese et al.1 
f2002b : r41 ISanderson et alJ42003[) ; I51 ILancaster et al.1 120051) . 

NOTE. — Uncertainties are 68% confidence. 
il Physical scale in kpc/i~' units corresponding to 1' at the cluster redshift. b SZE properties from 
AMiBA7 at 94GHz: cluster peak SZE flux (mJy) and angular size (') in FWHM measured from 
the cleaned image (Wuetal. 2008a). c Published X-ray properties: X-ray temperature (keV), X-ray 
core radius (kpc/r 1 ), and references. For A2142, Tx and 9 C are taken from Ref. [1], and Refs. [4,5], 
respectively. For A2390 a 10% error is assumed for (Tx, /3), for which no error estimate was presented 
in the original reference. 



TABLE 2 

Subaru weak lensing data and background galaxy sample 

Cluster Filters Seeing" n g b B/R c (D ds /D S ) A Zl , D e cr K f 

(arcsec) (arcmin -2 ) 

0.70 ±0.02 0.70+°^ 0.029 
0.88 ±0.04 0.95+° : ™ 0.021 
0.72 ±0.04 0.98+g:24 032 

0.72 ±0.04 1.00 + ii : ift 026 

a Seeing FWHM in the final co-added image in the redder band. b Surface number den- 
sity of blue+red galaxies. c Fraction of blue to red galaxies in the blue+red background 
sample. d Distance ratio averaged over the redshift distribution of the blue+red sample. 6 Ef- 
fective source redshift (see eq. 1221 ) corresponding to the mean depth {D ds /D s ) f RMS 
noise level in the reconstructed k map. 



A 1689 


Vi' 


0.88 


8.8 





A2142 


g'Rc 


0.55 


30.4 


2.1 


A2261 


VR C 


0.65 


13.8 


1.5 


A2390 


VR C 


0.70 


20.7 


2.1 



TABLE 3 

Summary of best-fit mass models from Subaru distortion data 



Cluster 


Tangen 

SIS 

<T V x 2 /dof 6»e 
fknu" 1 ) (") 


tial reduced shear, g+ 

NFW 

M vil Cyh X 2 / d °f #E 
(W i5 M Q /h) (") 


Lensing convergence, k 
NFW 

M vir c vlr X 2 /dof 9e 
(W ,s M Q /h) (") 


A1689 
A2142 
A2261 
A2390 


1403 ±41 11/9 47 ±3 
970 ± 27 39/8 25 ± 1 
1276 ±43 8.7/8 37 ±3 
1139±38 3.8/8 30±2 


7.3/8 47 + }5 
i 07+° :22 5 6 +(W 2 1/7 1 2+ 2 - 9 
1.35+°- 2 « 6.4+J' 7.7/7 20+jS 
0.9o3 : |I 6.9+f | 3.8/7 15+' 3 


1.05+Ojl 15.8+' 4 ,, 2 5.3/8 46+ 2 ^ 
I.06+™ 4.9+j-' 2 20/10 0.5+ 2 - 3 
1.26+<>- 2 ° 10.2+jj 9.8/8 37+ 25 , 
0.92+» : |5 7.3+f| 8.1/8 17+r5 



NOTE. — A flat prior of c V ir < 30 is assumed for the halo concentration of the NFW model. The Einstein radius 6e is calculated for a background source at 
z s = 1-5, corresponding roughly to the mean depth of blue+red background galaxies. 



TABLE 4 

Cluster mass models for gas mass fraction measurements 



Cluster 


NFW model 

Myir Cvir 
(lO 15 M0/r') 


SIS model 
(kms -1 ) 


''2500 ^500 100 r v j r 

(Mpc/r 1 ) (Mpc/r 1 ) (Mpc/r 1 ) (Mpc/r 1 ) 


A1689" 
A2142 
A2261 b 
A2390 


1 55+ 013 12 3+ 9 

1.07+a- 22 5.6+^ 

1 ?5+<' :i ' 1 1 1+ 22 
-8-16 -l»9 
Ogn+0 15 6 9+ 23 


1403 ±41 
970 ±27 
1276 ± 43 
1139 ±38 


0.57±0.01 1.16±0.02 1.70±0.04 2.13±0.05 
0.43±0.02 0.99±0.04 1.51±0.07 1.98±0.10 
0.52 ±0.02 1.06 ±0.04 1.56 ±0.06 1.94 ±0.07 
0.42 ±0.03 0.92 ±0.04 1.38 ±0.06 1.73 ±0.07 



a The NFW model is constrained by a joint fit to ACS strong lensing and Subaru distortion+magnification data, presented 
in Umetsu & Broadhurst (2008), but with our improved color selection of the red background sample for Subaru distortion 
measurements ( fl4-5.2l .k The NFW model is constrained by a joint fit to the inner Einstein-radius constraint and the outer 
Subaru k profile ( A4.5.21 . 
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TABLE 5 
AMlBA VISIBILITY ANALYSIS 



Cluster 


KS01 




isothermal /3( 


= 2/3) 
Y(3') 

(io- 10 ) 




yo 
(io- 4 ) 


F(3') 
(10- 10 ) 


(10- 4 ) 


A1689 


4.15 ±1.00 


2 5+°' 6 


4.31 ±1.10 


2.6+ 06 


A2142 
A2261 
A2390 


2.29 ±0.28 
3.00 ±0.84 
2.87 ±0.61 


3 5+°^ 
1 5+ - 5 

1 Q+0.6 


2.00 ±0.25 
4.25 ± 1.22 
3.40 ±0.72 


4 0+^ 
1 6+ 05 

7 1+0.8 



NOTE. — The effects of radio point source co ntamination in the ther- 
mal SZE have been corrected for (see Liu et al. 2008). The relativistic 
correction to the SZE is also taken into account. 



TABLE 6 

Cluster gas properties derived from the AMiBA/Subaru data 



Cluster 


KS01 + NFW 

M B as,2500 M gas ,500 Wgas,200 /gas.2500 /gas,500 /gas,200 

(lO 13 M /r 2 ) 


isothermal /3(= 2/3) + SIS 
Mgas.2500 Mgas,500 Msas.200 Aas.2500 /gas,500 /gas,200 

(io 13 M r 2 ) 


A 1689 
A2142 
A2261 
A2390 


4.4+j;' 8.8^| 11.5+1° 0.098+g;°| 0.115^| 0.11«$£ 

2.3 + ° A 3 7.2+J-5 11.2+2S 0.128+gMS 0.169+°-°*j 0.183«f« 

3.0+^ 63 + ' 2 \ 8.4+2] 0.087+°™ 0.103+™« 0.108+™ 
^ 3+0:7 61+2:4 q+4:o n 197+0:059 rii<;n+o:o75 n 1ft/1 +o:o84 

2.3_, s 6.1_, s 8.8_, 7 0.122_ nm1 0.153_ nnio 0.164_ nn „ 


3-SSS 7.8!2.0 H.8+3.0 . 100 +0.024 . 108 +O.026 O . m +0.027 

2.5+| : ? 5.4!j;5 8.1+| 3 0.097^|° 0.103+°|> 0.105+°|| 
2.2± 5.6^1 8.8-1! O.I25;'l!i 0.145^ 0.151$$ 



Note. — The derived gas fractions fggs scale wilh ihe Hubble parameter has /g^ oc /; (/; -0.7 adopted here). Confidence intervals are quoted al the 1<t (68%) level. Here we exclude ihe results from the isothermal model for A2142 which overpredicis 
/gas&t ull relevant radii (r > ''2500* com P' lr cd wilh the cosmic baryon fraction, = f2^,/f2 m = 0.171 ± 0.009. 
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FIG. 1 . — The quadrupole PSF anisotropy field for individual clusters as measured from stellar ellipticities before and after the PSF anisotropy correction. For 
each cluster field, the left panel shows the raw ellipticity field of stellar objects, and the right panel shows the residual ellipticity field after the PSF anisotropy 
correction. The orientation of the sticks indicates the position angle of the major axis of stellar ellipticity, whereas the length is proportional to the modulus of 
stellar ellipticity. A stick with the length of 5% ellipticity is indicated in the top right of the right panel. 
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FIG. 2. — Stellar ellipticity distributions before and after the PSF anisotropy correction for individual clusters. For each cluster field, the left panel shows the 
raw ellipticity components (e* , e|) of stellar objects, and the right panel shows the residual ellipticity components (5e* , <SeJ ) after the PSF anisotropy correction. 
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FIG. 3. — Top panels: Mean shape distortions (g+.gx) averaged over the entire cluster region (1' < < 18') for the four clusters done separately for the 
blue and red samples, in order to establish the boundaries of the color distribution free of cluster members. Bottom panels: Respective numbers of galaxies as a 
function of color-limit in the red (right) and the blue (blue) samples. 
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FIG. 4. — Mass maps of the central 22' X 22' of four AMiBA/Subaru clusters reconstructed from Subaru weak lensing data, with the gravitational shear field 
of background galaxies overlaid; 10% ellipticity is indicated top right, and the resolution characterized by Gaussian FWHM is shown bottom right. Also overlaid 
are contours of the SZE flux densities at 94 GHz, observed with the 7-element AMiBA, given in units of \a reconstruction error. The resolution of AMiBA, given 
in Gaussian FWHM, is 6' For all four clusters the distribution of the SZE signal is well correlated with the projected mass distribution, indicating that the hot gas 
in the clusters traces well the underlying gravitational potential dominated by unseen dark matter. The dark blue regions in the mass map of A2142 are outside 
the Subaru observations. 




FIG. 5. — Azimufhally-averaged radial profiles of the tangential reduced shear g+ (upper panels) for the four clusters based on the combined red and blue 
background samples. The solid and dashed curves show the best-fitting NFW and SIS profiles for each cluster. Shown below is the 45° rotated (x) component, 

gx- 
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FIG. 6. — Model-independent radial profiles of the lensing convergence k(8) = £ m (t9)/X C rit for the four clusters derived from a variant of the non-linear 
aperture-mass densitometry. For each cluster, the best-fitting NFW model for the re profile is shown with a solid line. The dashed curve shows the best-fitting 
NFW model for the g+ profile in Figure[5] 
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FIG. 7. — Subaru V +R C pseudo-color image of the central 6.7' X 6.7' (2K X 2K pixels) region of the cluster A2261 at zd = 0.226. Overlaid is the tangential 
critical curve predicted for a background source at z, s ~ 1.5 based on strong lensing modeling of multiply-lensed images and tangential arcs registered in deep 
Subaru VR C and CFHTAVIRCam JHK$ images. The effective radius of the tangential critical curve defines the Einstein radius, Be rs 40' at z s ~ 1.5, 
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FIG. 8. — Cluster ellipticity and orientation profiles on mass and ICM structure as a function of aperture radius t9 ap , measured from the Subaru weak lensing 
and AMiBA SZE maps shown in Figure[4] For each cluster, the top panel shows the halo ellipticity profile |e halo |(c9ap), and the bottom panel shows the orientation 
profile 20((9 ap ), where </> hal ° represents the position angle of the major axis as measured from weighted quadrupole shape moments. 
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FIG. 9. — Gas mass fraction profiles (/gas(< r)) = (M gas (< r)/M tot (< r)) averaged over the sample of four hot (T x > 8keV) clusters (A1689, A2142, A2261, 
A2390) obtained from joint AMiBA SZE and Subaru weak lensing observations, shown for the NFW-consistent Komatsu & Seljak 2001 model (black) and the 
isothermal /3 model with fi = 2/ 3 (bl ue), along with published results (square, triangle, and circle) from other X-ray and SZE observations. The isothermal results 
exclude the cluster A2142 (see j|6.4l . For each model, the cross-hatched region represents lcr uncertainties for the weighted mean at each radius point, including 
both the statistical measurement uncertainties and cluster-to-cluster variance. The black horizontal bar shows the constraints on the cosmic baryon fraction from 
the WMAP 5 -year data. 
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APPENDIX 

ONE-DIMENSIONAL MASS RECONSTRUCTION FROM DISTORTION DATA 

Following the method developed bv lUmetsu & Broadh urst ( 2008), we derive an expression for the discrete convergence profile 
using a non-linear extension of weak lensing aperture densitometry. 

Non-Linear Aperture Mass Densitometry 

For a she ar-based estimation of the cluster mass profile we use a variant of weak lensing aperture densitometry, or the so-called 
(-statistic dFahlmanet al.lll994t IClowe et alfeOOOl) of the form: 



C c (e) = 2 / d\n9' 1+ {9') 
Je 



d\n0' 1+ (6') 



+ l-(0inn/0out) 2 A»„ 

= R(9)-Me mn <d<e oat ) 1 (Ai) 

where n{9) is the azimuthal average of the convergence field k(0) at radius 9, R{9) is the average convergence interior to radius 
9, #; nn and out are the inner and outer radii of the annular background region in which the mean background contribution, 
Rb = R(9 mn < d < fl ou t), is defined; i+ (9) = R(9)- k(9) is an azimuthal average of the tangential component of the gravitational 
shear at radius 9 (Fahlm an et al.l[l9 94), which is observable in the weak lensing limit: "f+(9) « (g+(9)). This cumulative mass 
estimator subtracts from the mean convergence R(8) a constant Rb for all apertures 9 in the measurements, thus removing any DC 
component in the control region 9 = [9i nn , 9 0Ut ]. Note that the Rb is a non-observable free parameter. This degree of freedom can 
be used to fix the outer boundary condition, and hence to derive a convergence profile k(9). 

In the non-linear regime 7+ (fl) is not a direct observable . Therefore, non-linear corrections need to be taken into account in the 
mass reconstruction process ( Ume tsu & Broa dhurst 2008). In the subcritical regime (i.e., outside the critical curves), 7+(#) can 
be expressed in terms of the averag ed tangential reduced shear as (g+(9)) ~ 7+(fl )/[l - k(9)] assuming a quasi-circular symmetry 
in the projected mass distribution (Broadhurst et al. 2005b; fOmetsu et ai1l2007l) . This non-linear equation (IA1 b for ( c (9) can be 
solved by an iterative procedure: Since the weak lensing limit {k, \ \g\ <C 1) holds in the background region ft nn < 9 < 9 max , 
we have the following iterative equation for Cc(9)- 



C ( c k+l) (9)^2 / dln9'{g + (9))[l-K«\9)] 
Je 



d\n9'(g + (9')), (A2) 

i ~ \ Vmn/ Vout) J iun 

where Cc* +1) represents the aperture densitometry in the {k+ l)th step of the iteration (k = 0, 1,2, ...,A^ ter ); the is calculated 
from Cc* +I) using equation ( lAlOt . This iteration is preformed by starting with k (0) = for all radial bins, and repeated until 
convergence is reached at all radial bins. For a fractional tolerance of 1 x 10~ 5 , this iteration procedure converges within A^ ter ~ 10 
iterations. We compute errors for ( c and k with the linear approximation. 

Discretized Estimator for the Lensing Convergence 
In the continuous limit, the averaged convergence R(9) and the convergence k(9) are related by 



2 f e 



dln9'9' 2 K(9'), (A3) 



For a given set of annular radii 9 m (m= 1,2, ...,N), discretized estimators can be written in the following way: 

2 m_1 

K m = R{9 m )=— ^Aln0,0>(0O, (A5) 
m i=i 



where 



ni = K{9i) = a 2 K M -a x Ki (/ = l,2,...,jV-l), (A6) 



1 2Aln6», \9,J ' " z 2Aln6>, V 9, 



with A In 9 1 = (ft +I - 9/)/ 9/ and 9/ being the area-weighted center of the /th annulus defined by 9[ and 6>/ +1 ; in the continuous limit, 
we have 



9, = 2 d9'9' 1 /{9] +l -9 1 l ) 
Je, 
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The technique of the aperture densitometry allows us to measure the azimuthally averaged convergence R(0) up to an additive 
constant Rb, corresponding to the mean convergence in the outer background annulus with inner and outer radii of 8 mn and 8 out , 
respectively: 

S(0) = + (A9) 
Substituting equation (IA9t into equation (IA6l l yields the desired expression as 

<e l ) = a l 2 ae M )-a\U0i) + (a l 2-a\)R b . (A10) 
Finally, the error covariance matrix of K; is expressed as 

-c^a l 2 Cl M -c/ 2 a[Q +hl) (All) 

where c£ = (S^Ci) i s the bin-to-bin error covariance matrix of the aperture densitometry measurements which is calculated by 
propagating the rms errors ct+(6*/) for the tangential shear measurement. 



